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Executive Summary 


This technical report (quarterly) details the work for Office of Naval Research (ONR) by Tennessee Tech. The goal 
of this project—jointly funded by ONR, NSF, and ARO—is to build a general purpose testbed with time reversal 
capability at the transmitter side. The envisioned application is for UWB sensors and tactical communications in 
RF harsh environments where multipath is rich and can be exploited through the use of time reversal. The report 
summarizes the results for each of two major tasks. Specifically, the project can be broken into 


• Task l-Theoretical Research 

• Task 2-Experimental Testbed 

[n the part of theoretical research, the central result is to enable nonTading transmission through the use of time re¬ 
versal Next, the optimum pulse shaping is also considered, to extend the range. Finally, the chirp pulse waveform— 
practical and useful—is considered. This chirp UWB system, combined with lime reversal, will be promising for 
long range sensor networking. 

In the experimental testbed part, progress in the second generation of the testbed—with time reversal—is described. 
At this moment of writing, the first generation of the testbed—without time reversal—is working at the Pi's Lab. 
over the air with a bandwidth of more than 500 MHz. The second generation will be built upon the first generation. 
The primary modification is the arbitrary modulation waveform generation at the transmitter. 
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Chapter ! 


Introduction 


Civilian and naval vessels have long been potential targets for criminal and terrorism activities. There are concerns 
that terrorists can ship various types of weapons of mass destruction (WMD) to international ports using commercial 
ships and their cargoes. One of the major steps in preparedness for such danger is reliable wireless communications 
both between the boarding party as well as effective ship-lo-ship and ship-to-shore communications once the first 
sign of threat is detected. This need is made possible, only recently [I], with the advent of revolutionary ultra- 
wideband (UWB) technology [2, 3, 4, 5, 6 3 7, 8, 9]* This proposed research is motivated for RF challenged envi¬ 
ronments such as hospital, coal-mine, intra-ship, intra-vehicle, intra-engine, manufacturing plants, assembly lines, 
nuclear plants, body area network sensors surrounded by vehicles and tanks, etc. 

In the part of theoretical research, the central result of this report is to enable non-fading transmission through the 
use of time reversal. Next, the optimum pulse shaping is also considered, to extend the range. Finally, the chirp 
pulse waveform—practical and useful—is considered. This chirp UWB system, combined with time reversal, will 
be promising for long range sensor networking. 

In the experimental testbed part, progress in the second generation of the testbed’—with time reversal—is described. 
At this moment of writing, the first generation of the testbed—without time reversal-—is working at (he Pi's Lab, 
over the air with a bandwidth of more than 500 MHz. The second generation will be built upon the first generation. 
The primary modification is the arbitrary modulation waveform generation at the transmitter. It is expected that the 
receiver part of the second generation mainly reuses that of the first generation. 

A key understanding is achieved so far: to further extend the range, the design of modulation waveforms must 
be carefully selected, according to some criteria. A good criterion is to maximize the signal to noise ratio (SNR) 
measured at the point after the matched filter. The frequency selectivity caused by the huge bandwidth (or pulse 
waveform distortion in the absence of multipath) will lead to the unique problem of pulse waveform optimization. 

Another key understanding is that the effective channel enabled by time reversal can be treated ‘‘Non-Fading”—a 
dream for wireless engineers from half an century. The system design will greatly simplified under this non-fading 
framework. Although this theoretical breakthrough ts reached, the design principle and system parameters to exploit 
this new phenomenon are far from clear. 

One lesson learned from this research is that time reversal is basic, for a dense multipath environment—a continuous- 
time impulse response of extremely long delay spread. 
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CHAPTER L INTRODUCTION 


The key UWB system design consists of two parts: (1) time reversal to deal with dense multipath; (2) the selection 
of UWB modulation pulse waveform—software defined or cognitive. From a information-theoretical viewpoint, the 
two parts as a whole form so-called “pre-coding”. 
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Chapter 2 


Non-Fading Transmission 


We propose a time reversal system paradigm for localization, communication, and networking. UWB transmission 
may lead to (almost deterministic) no-fading communication—a dream for wireless industry over half a century—to 
enable some applications, for which narrowband techniques may fail. 

The transceiver structure has been the central topic in the UWB community. The Pi’s research roughly follows 
this trend. First, the RAKE structure [10, II, 12, 13, 14, 15, 16, 17] is used, and found too costly—after a long 
painful search, due to the large number of paths (fingers) and timing requirement [6, 7, 8, 9], In particular, pulse 
waveform distortion will cause a design problem, first recognized (one decade ago) by the PI based on his channel 
model [10, II, 12, 13, 14, 15, 16, 17], which was accepted into IEEE standard [18], The orthogonal frequency 
division multiplexing (OFDM) system [19,20] is slightly more complicated than the RAKE system. Secondly, the 
transmitted reference (TR) structure [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31* 32] is used to reduce the timing 
requirement, and avoid the need of channel estimation. Thirdly, for sensor networks of low-cost and low power, 
non-coherenl energy detection (analog frontend) [33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43] is preferred to the TR, 
although both systems are adopted in IEEE 802J5.4a. Inter-symbol-interference (ISI) limits the achievable data 
rates of the TR and energy-delection systems. 

Fourth, for data rate higher than the TR and energy-detector, time reversal at the transmitter—perhaps viewed as 
continuous-time, channel-matching filter, or generalized RAKE [17] at the transmitter—is combined with transmit¬ 
ted reference structure [44, 45, 46, 47], and non-coherenl energy detection [48, 49, 50]. It is even extended to the 
UWB multiple input and multiple output (MIMO) system [8, 52, 53, 54, 55, 56, 57, 58]. Time reversal has been 
used in UWB radio [59, 60, 61,62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73]. It is used extensively in acoustics 
by Fink’s team [74, 75,76,77, 78, 79, 80, 81, 82, 83, 84, 85,86], and later to underwater acoustic communications 
mainly by Kuperman’s team [87, 88, 88, 89, 90, 91,92, 93, 94, 95, 96, 97, 98, 99, 100, 101, J02* 103] to deal with 
rich multipath phenomenon which is similar to UWB. It has been first used in wire-line communication by Bell 
Labs [104, 105] and IBM [106]. In spirit, time reversal is similar to the classical phase conjugation concept in 
optics [107, 108, 109], and retro-directive array in microwaves [110, Ml, 112, 113, 62, 114, 115, 116, 117, 118, 
119, 120, 121, 122, 123, 124, 125, 126]. The objective of ail the above schemes is to exploit the channel itself as 
part of channel-matched transceiver. This principle has the advantage of real-time, since the operation is finished 
in a continuous-time, analog form, in contrast to the algorithm convergence of the discrete-time form used by such 
systems as equalizers [127, 128, 129, 130, 131, 132, 133]. A key observation [53, 52, 134, 135, 136) has been 
recognized that, if lime reversal is always used at the transmitter side, no fading needs be explicitly addressed in a 
UWB impulse radio. This effect derives from the (multipath and frequency) diversity gain offered by the huge signal 
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CHAPTER 2. NON-FADING TRANSMISSION 


bandwidth. The proposed research is based on this key cffccl—to exploit multipath and signal bandwidth (or pulse 
duration). 

With time reversal—of the channel impulse response (CIR) illustrated in Fig. 4.3—at the transmitter side, the 
effective channel between a pair of a modulation waveform and a demodulation waveform is the autocorrelation o* 
the CIR—illustrated in Fig. 2.1(c), This new channel is much more stationary than its old counterpart [134, 135]. 
An analogy is as follows. If the CIR is a realization of a random process—the first order statistics is changing 
rapidly from one realization to another realization, the second order statistics is much more reliable and stable. Its 
time-varying properly highly relies on the number of multipath which* in return, is determined by the pulse signal 
bandwidth, as a function of RF environments. The significance of this proposal is to convert the conventional CIR 
problem into a new, autocorrelation problem, and systematically establish a theoretical framework based on this 
principle. This understanding is in parallel with the historical development of the white Gaussian noise theory in 
early 1950s [137, 138* 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151,152]: what really matters is 
not the CIR’s first order statistics, but its second order statistics! 

Once the key conceptual breakthrough is made, the formulation of this idea in a mathematical framework is rela¬ 
tively straightforward: The system model consists of a lincar-timerinvariant (LT1), dispersive channel with additive 
(Gaussian or non-Gaussian) noise and interference [20]. We have converted a time-varying, dispersive (wireless) 
channel into a time-invariant, dispersive channel—that is much easier to handle. An exhaustive research has been 
done regarding this classical system model [153, 154, 19, 155, 156, 133, 157, 158, 159, 160], since it is the model 
used for wire-line ! 153, 156, 133, 157, 19]—DSL and VDSL [19]—for which the Bell systems have the deep-packet 
resources [127, 128, 129, 130, 131, 132, 133], The well-understood model is an excellent tool for the proposed 
research, and paves the way for our next engineering breakthrough. As one example, the modal modulation used 
by Holsinger (1964) [153] is a starting point for our new adventure. Time reversal is used at the transmitter. The 
so-called time reversal modal modulation is proposed here. Our system model consists of an effective LTI, dispersive 
channel 1 with additive (Gaussian or non-Gaussian) noise and interference. The simplified mathematical framework 
will lead to a better design, since the dimensions of bandwidth and multipath can be focused on, through the effective 
LTI channel. The unified theory links up all the fundamental system parameters such as bandwidth, multipath, and 
space, and box them into an effective CIR. In practice, all we only need to do is to replace the LTI model of the clas¬ 
sical theory with our effective LTI model: it is that simple! Now, the open question is to spell out all the necessary 
conditions—such as pulse duration and waveform as a function of RF environments—for this new framework to be 
valid. 

Built upon the above new understanding, theory and prototype testbed are proposed to lest and validate this new 
system design paradigm. A 5 GHz UWB radio testbed has a special meaning in the context of the new paradigm, 
since, for indoor channels, 2 GHz has been found to be the critical bandwidth, beyond which the effective channel 
can be considered as “deterministic (non-fading) in practice, provided receivers exploiting the full channel energy 
are employed” quoted from [161 f. Mixed signal processing is the bottleneck* To overcome this problem, a new 
patented design [162, 163] from radar community can be used. 

Similar to short pulse radar [164, 165, 166, 167, 168, 169], waveform optimization for detection is a basic problem. 
For UWB MIMO, we have derived a formulation based on the eigenfunctions of the integral equation. This problem 
may be related to the multiuser MIMO problem that can be effectively solved via a convex optimization [170, 
171, 172, 173, 174, 175 L Bear in mind the fact that, since a practical UWB system works in a very low SNR 
range, optimization for detection (maximizing SNR) is the primary problem—not the one to maximize the mutual 

'if the CIR of the LTI, dispersive channel is h(r), then the CIR of the effective LTI channel is h( — r) * h(r) where is the linear 
convolution, 

: in the long process of studying brae reversal, the Pi’s team recognized this result, independently of f 161). 







2.1. SIGNIFICANCE 


9 


information. 


2.1 Significance 


Time-varying fading is the hallmark of narrowband wireless communications. It seems safe to say that almost all the 
major progress made in wireless industry, such as MIMO and OFDM, has been motivated to exploit fading through 
(multipath and spatial) diversity. Fading caused by multipath, however, may be negligible for UWB communication 
using a sufficiently large signal bandwidth, e g. 2 GHz for indoor [161]. The proposed research is motivated to 
explore the dimensions of bandwidth and multipath in real-life environments—an issue unique to UWB impulsive 
radio. The number of paths is roughly proportional to the bandwidth, when it is below 2 GHz for indoors. In 
particular, an extremely large number of (time-invariant) resolvable multiple paths—treated as continuous-time, 
dense waveform over a duration of 800 ns, see Fig. 4.3—will be encountered in real-life applications [l, 176, 177, 
178,179], and will be investigated as an example. When fading can be ignored in practice, the mathematical tools are 
simplified to a level accessible by engineers with only undergraduate education [20, 180]. The body of knowledge 
based on fading communication mathematics [181, 182]—generally very difficult for engineers and students, even 
at a graduate level—can be circumvented, perhaps for the first time, in the wireless communication curriculum. The 
new mathematical framework will pave the way for a novel system paradigm and a new understanding of the mystic 
multipath. 


2.2 UWB Short Pulse Propagation 

The non-fading (deterministic) treatment of a wireless channel is based on the experimental study of the UWB 
propagation channel. There is a class of extremely rich multipath environments. This class of channel has been 
recently investigated by the Pi’s team [53, 52, 134, 56, 57, 58, 135]. It is found that the number of paths is much 
larger than that of the office environment—compare Fig. 4.3 and Fig. 2.1(d). Communications and networking in 
confined metal environments [134, 135, 57] in this class—Fig.2.1(a)—is motivated for applications such as intra¬ 
ship [1], intra-vehicle [176], intra-engine [52], manufacturing plants, assembly lines, nuclear plants, body area 
network sensors surrounded by vehicles and tanks, etc. 


2.2.1 Related Work 

One difficult problem is communications over such RF challenged environments as a confined metal box. Narrow 
band wireless technologies have proved ineffective in these environments, due to resonance caused by the metal 
walls [177, 178, 176,179]. Chirp-like signals are observed in a jet engine [178]. Only recently [ I ], this problem has 
the promise to be solved with the advent of UWB technology. Lawrence Livermore National Laboratory (LLNL) 
just demonstrated UWB transmission in a ship in 2006 [ I ], The PFs team demonstrated this in an engine [52J, and 
an metal cavity [134, 135, 57]. Strict physics for the latter case is given by L. B. Felsen [177] who introduced the 
UWB problem to the PI in 1993, and W. Geyi [179] who was introduced to this problem by the PI in 2007 IEEE 
APS conference. In general, this problem is open, especially in the context of non-fading transmission. The mobility 
of antennas needs to be considered. 





10 


CHAPTER 2. NON-FADING TRANSMISSION 



Figure 2.1: UW13 short pulse propagates in RF harsh environments, (a) Rectangular metal cavity: (b) OR in a rectangular 
metal cavity. The pulse has its spectrum from 3 GHz to 10 GHz • (c) autocorrelation of CIR in (b); (d) CIR in an office 
environment. The experimental setup is identical to that of (b). 


2*2.2 Research Plan 

The central task of the proposed research in this context is to validate the system design assumption, and parameterize 
the effective channel h(— r) * fr(r). When the testbed is available, the field test will be carried out, to demonstrate 
the system principle in different RF harsh environments. A lot of new experiments are needed, in the process of the 
proposed research. Fortunately, the PI, a veteran in channel modeling, has the expertise and the required equipment 
to perform these tasks. 


Preliminary Work 

The PI has accumulated more than one decade’s experience in this channel modeling area. His Lab has two the- 
state-of-the-art channel sounding systems for UWB communication. The time domain system Fig, 2.3(b) is based 
on digital sampling oscilloscope (DSO), Textronix TDS 7000E3, The time domain sounder has the capability to 
handle the MIMO channel sounding—not using virtual array approach. The frequency domain system uses a vector 
network analyzer (VNA), Agilent N5230A (300kHz-13.5GHz). These sounders helped to produce two PhD degrees, 
five MS degrees, and five undergraduates used them through NSF s REU (Research for Undergraduates) program. 

Fig. 2.1 illustrates a new phenomenon—an extremely long train of multipath echoes. Fig. 4,3 contains much more 
paths than that of Fig. 2. I (d). The measurement is performed using frequency domain technique to analyze the 
characteristics of the UWB channel in the confined metal environment. The rectangular confined metal cavity in 
Fig. 2J(a) measures 16 feet by 8 feet by 8 feet. The materials of the walls are aluminum. The sounding is carried 
out by sweeping a set of narrowband sinusoid signals (tones) through a wide frequency band. The data of channel 
frequency response is collected using the VNA based channel sounder in the PTs lab. I he VNA sweeps from 3 GHz 
to 10 GHz using 7001 points with frequency step of 1 MHz. The power of each tone is 10 dBm. 


2.3 Time Reversal Modal Modulation 

2.3.1 Related Work 

When time reversal is used, the effective channel is regarded as deterministic (no-fading) channel. The proposed 
time reversal modal modulation is different from OFDM in several ways. First, OFDM is a special limit of this 
modulation, when the pulse duration becomes infinite. Most UWB pulses are short in duration, similarly to UWB 
























2.3. TIME REVERSAL MODAL MODULATION 




Figure 2.2: Time reversal modal modulation is proposed to simplify the transceiver complexity. The channel matched filler 
h(T h -t) is moved to the transmitter side. The self-reproducing can be delined. 

radar. Secondly, OFDM can be viewed as multiple tones (multiple carriers), implying each tone assumes a single 
frequency value. This is not true for the proposed modulation. The effective channel can be diagonalized by its 
eigenfunctions (eigen-waveforms)—changing with RF environments, similar to the radar problem [165, 166, 167, 
168, 169], These eigen-waveforms arc used for modulation, and may be wideband (say 100 MHz). Thirdly, the 
pros and cons of this modulation relative to OFDM are unclear, in the context of system implementation* 

23.2 Research Plan 

The proposed system formulation for optimum modulation waveform is based on two assumptions: 


* The channel is linear, time-invariant (LTI) for every symbol transmission. 

• The spectra of noise and interference signals are available at both the transmitter and receiver* 

The first assumption has been established recently by the Pi's group, through experiments [53, 52, 134, 56, 57, 58, 
135], in the framework of time reversal. Knopp’s team reaches the same result [161 ], from a channel modeling view. 
The second assumption is the practice in current industry. 

The mathematical formulation for the so-called modal modulation traces back to the well known work oflfolsinger 
(1964) at MIT [153, 19], originally motivated for telephone line—widely regarded as the first work for OFDM— 
regarding optimum modulation for a deterministic, dispersive, LTI channel with memory. It is found that Of DM is 
an asymptotic case of Holsingeris modal modulation, when the duration of the modulation pulse goes infinite. f : or 
a nanowband technique, the modulation pulse is a sinusoidal pulse of infinite duration* Since modal modulation 
is optimum—to be made more clear later, OFDM is accepted widely without question. With the advent of UWB 
devices, we must reexamine the practice for narrowband techniques. With the aid of the effective channel introduced 
in Section 2, the mathematical framework of Holsinger’s is sufficient: Our task is to replace his LTI channel mode] 
with the effective channel model—comparing Figs, 4,6 and 4*7, 

It is found recently by the PFs team that modal modulation is not equivalent to OFDM for UWB modulation pulses 
of short finite duration. It appears that this finding requires a paradigm shift in that UWB transmission does not have 
to use a sinusoidal pulse! 

Single Input and Single Output (SISO) Let us fix some notations for more precise discussions about Ho I singer's 
modal modulation, illustrated in Fig* 4*6, for colored noise and an arbitrary, finite observation interval. Ifx(*) is 
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the input to the LTI channel—based on assumption I—with channel pulse response /»(/). The symbol duration of 
x(t) is [0,7 s], the arbitrary, finite waveform duration of iff) is [0,7'j—which implies that x(t) is equal to zero 
if t is between 7 and T$. The output waveform r (?) ~ J f| x (r)h (t — r) dr , 0 < t < /i, for an arbitrary, finite 
observation interval [0,7]], A colored noise with a spectral density N(f) —based on assumption 2—defines a kernel 
function K\ (f, s) by 

R n (t — a)K\(a 7 s)da — 6(t — s), 0 < t,s < I\ (2.1) 

where /^(r) — N(f)eJ uir 4f defines a function K(t, s) by 

K(t r s) = f j h(a — i)K\(cr 7 s}h(p — t)dadp 0 < t f s < T (2.2) 

JQ JO 

and define a possibly infinite-size set of eigenfunctions, through a integral equation [183, 184] 



An • <fn(t) ~ C K(t t T)Vn(T)dT, 0 < t < T (23) 

Jo 

That is a function <^(£) that when convolved with K(t 7 s) over the interval [D t Tj reproduces itself, scaled by 
a constant A n , called eigenvalue. The Pn{t) n = 1,2,..., oo form a set of complete, orthonormal basis func¬ 
tions, so x(t) —an arbitrary modulation waveform of finite duration T—can be expanded in terms of x (/) - 

52 x n*Pn (£), 0 < t < T and x n — / (( x (7) <p n (t)dt, n — 1,2,**. ,oc, so the energy of the transmitted sig- 
n-0 

nal waveform, E t — J { ( [x (t]| 2 d£, can be expressed as tit ~ 52 a: 2 * When while Gaussian noise is assumed, a 

n—1 

00 

simple matched filter can be used: Hie output signal to noise ratio (SNR) is 2E/!'i\, where E — 52 Aand 

Tl“l 

Nq/2 is the double sided noise spectral density. So the probability of error for the detector is given by 





(2,4) 


In order to minimize / J e , SNR should be maximized. Since, by convention, Ai > A 2 > A-$ > * * ■, if the energy of 
is fixed, maximizing E means a:(£) — —the optimal modulation waveform. So, the maximum SNR 

2A x^ 

can be expressed as SNRmax = — ^ 1 with the maximum eigenvalue A^ — A] for our above notation. This is 
ideal for range extended UWB sensors—the primary goal of this project. 


From another (capacity) point of view, the modal modulation system have many parallel orthogonal channels, and 
each channel corresponds to one eigenfunction r jh(0- We can use the classical water filling scheme to calculation 
coefficients x n n — 1 , 2, ... ,oo, based on the energy constraint of transmitted signal x(r). 

A new time reversal modal modulation, illustrated in Fig. 4.7, is proposed to simplify the transceiver complexity. 
The h(l\ — /), time reversal filter, is moved to the transmitter side. Without loss of generality, wc use white Gaussian 
noise as an example to illustrate the concept. Similarly to the modal modulation, the self-reproducing eigenfunctions 
4> n (t) with eigenvalues k„ illustrated in Pig. 4.7 can be defined as 

tin ■ <t>n(t) = [ Rhh{t - r)^ n (r)dr, 0<t<T (2.5) 

Vo 


where Rhh(t) ™ A(£) * A{—7). There are four proposed tasks in the SfSO system: 
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Figure 23: UWB MIMO System, (a) MIMO System: (b) Channel Sounding, 


• Concept proof The no-fading mathematic framework needs to be proved, and made more precise, by spelling 
out system parameters such as bandwidth and waveform. Real-world experimental data will be used for this 
purpose, 

* Optimization of modulation waveforms over different RF environments This requires modulation wave¬ 
form optimization, as a function of waveform duration T, defined in Eq. (4.5). 

m Interference spectrum estimation and its impact Cognitive radio must have the function of “detect and 
avoid” other narrowband radios such as Wi-Fi and WiMax, This is a challenging problem, even for a narrow- 
band radio. One simple method is to estimate the background interference over the ultra-wide bandwidth—say 
5 GHz—and optimize the waveform, based on the key equation (4.5). 

Multiple input and Multiple Output (MIMO) Let us consider the system block diagram of the time-reversed 
impulse MIMO in Fig. 23. Let us denote the channel impulse response (CIR) relating the m-th element at 

the transmitter to the n-th dement at the receiver. The set of impulse responses is called the propagation operator 
If pulse waveform vector a(f) = [a\ (t), 03 (f), ...*<»a*(*)]* of finite duration /' is sent from the transmitter, then the 
vector b(t) = [fej * s received at the receiver, for an arbitrary, finite observation interval [ 0 , £]. 

The superscript H f” represents the conjugate transpose (Hermitian) of a vector or matrix. It follows that 

b(() = H(0 * *(t), 0 < t < T 2 (2.6) 

where H(f) is the matrix of N x M with elements of h nm (t ). The notation of represents element-by-element 
convolution. On the other hand, due to the spatial reciprocity, if one sends a signal n — 1 ,N from die receive 
array, the signal / m (f), m — 1 , M is obtained in the transmit array. Defining column vectors e and f as a and b, 
respectively, it follows that 

f(t) = H f (f) * e(f), 0 < t < r 2 (2.7) 

As a consequence, H (t.) permits one to calculate the forward propagation of impulses from the transmit array to the 
receive array, while Ht(t) the backward propagation from the receive array to the transmit array. If the time reversed 
matrix t) is used as part of the precoding matrix—in analogy with Fig. 4.7, the equivalent matrix channel is 

Rhh( 0 = H(t) * H+(-i) (2.8) 

Eq. (2.8) will contain the temporal-spatial focusing—unique to UWB pulse signals. One objective UWB MIMO 
is to enable non-fading transmission. Cooperative distributed antennas can be used. MIMO multiuser networking is 
our future research focus, in the framework of time reversal. 

in analogy with (4.14), the self-reproducing eigenfunction vector #,(f) can be defined as 

T* 

Kn - 0 n {t) ~ f R H h(* - r)e n (r)dr t 0 <t< T ( 2 . 9 ) 

Jo 

Once Eq. (2.9) is solved, the rest of signal processing is similar to the SISO system above. An alternative approach to 
solve for optimum MIMO waveform may use the recent result of information theory, based on convex optimization 
[170, 171, 172, 173, 174, 175]. There are three proposed tasks in the MIMO system: 
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Figure 2*4: Comparison of modulation waveforms for optimum modal modulation for different environments, (a) wedge 
(street comer); (b) rectangular plate (building); (c) waveform comparisons: A passband channel has an ideal passband filter as 
its frequency response; (d) optimum modal modulation waveform calculated using OR of Fig* 4.3 


• Concept proof Similar to the STSO system, experimental and theoretical work is required to extend the frame¬ 
work to the MIMO system* The focus is on the validity of the framework, by considering multiple antennas* 
The PFs lab has two channel sounding systems (in both time—in Fig* 2.3(b)—and frequency domain) for this 
purpose* 

• Spatio-temporal optimization Optimization of waveform is required, by taking into account the space di¬ 


mension, through Eq. (2.9)* 

• Spectrum estimation Sensing the spectrum is a critical part of cognitive radio. Algorithms for spectrum 


estimation need be investigated. 


2,3*3 Preliminary Work 

Compared with the benchmark of the time reversal chirp system, there is a big room for optimization (Fable ??)— 
as a function of waveform duration—by using time reversal modal modulation system. This has been justified 
from two different metrics. Different propagation environments define different optimum modulation waveforms, as 
illustrated in Fig* 2*4. 

From Detection Point of View The SNR after matched filter in the receiver side is the metric to determine the system 
performance. 

From Data Rate Point of View In the time reversal modal modulation system, there are many parallel orthogonal 
channels, and each channel corresponds to the orthonormal eigenfunction (f) and eigenvalue from Eq. 4.14* 
One waveform design example of the 4>n(t) is given in Fig* 2*4. We assume K| > > /s& > - ■ ■ * Meanwhile, 

the bits transmitted in these different channels are assumed to be totally independent* If L-level pulse amplitude 
modulation (PAM) is considered, there are L equally likely amplitude values. In time reversal system, the transmitter 
waveform filter is G T (t), while die time reversal filter is /*{-<)* In the time reversal chirp system, <3r(f) is chirp 
waveform. When BER P* is set to 10 4 and the same PAM constellation is considered, e.g*, for L — 8, the data rate 
gap (in dB) between time reversal modal modulation system and time reversal chirp system is shown to 7.5, ILL 
13*9, for the waveform duration of T (in ns) equal to I, 10, 100. 
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Chapter 3 


Optimum Transmission Waveform in 
Canonical Channels 


'flic purpose of this chapter is to introduce readers the time domain modeling for the UWB signals. For wide¬ 
band signals such as UWB* by principle, time domain frame work appears to be a better frame work to study its 
propagation. The time domain impulse responses for a series of scattering environments have been derived for the 
first time. Specifically, several electromagnetic communication fundamental channels, including perfect conduct 
electromagnetic wedge, large rectangular plates, and parallel planes waveguide, are studied. The receiving signal, 
calculated through convolution in the time domain impulse response, can provide waveform information, as well 
as the amplitude information (Path loss). The pulse waveforms passing through these channels are compared with 
their original waveforms. The comparison show that for wideband signal, pulse could be extremely distorted by 
diffraction mechanism. The results converted from frequency domain by applying a inverse fast Fourier Transform 
(IFFT) serves as sanity checks for our time domain derivation. Given the time domain channel impulse response, we 
also study the optimum transmit waveform such as the signal to noise ratio (SNR) at the receiver side is maximized, 
under the constraint of fixed transmitted power. 


3.1 Introduction 


The purpose of this chapter is to derive the channel impulse responses for a scries of environments and analyze 
the corresponding optimum transmit pulse waveforms. The UWB system are generally based on the transmission 
of many consecutive pulses. Therefore, by principle, time domain framework appears to be a belter framework to 
study the radiation, propagation and reception of UWB signals. Moreover, comparing with traditional frequency 
domain framework, lime domain framework provides more physical insight, ‘"or example, except predicting path 
loss, time domain convolution with TD-UTD can also predict the receiving pulse waveform and show the possible 
pulse distortion caused by diffraction, which has been proved to be a significant issue for wideband signals such as 
UWB [l]-[3]. 

The essential benefits of the time domain framework motivate this work. However, the traditional time domain 
methods such as Finite Difference Time Domain (FDTD) and time domain integral equation become intractable 
when the width of the incident pulse is very narrow compared with the geometric dimensions of the scattering 
objects, which is the case when we study the propagation of UWB signals. In this chapter, we will try to attack the 
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propagation problems for wideband signals directly in the time domain, with a closed form solution. Wc derived 
the time domain impulse response for a series of scattering environments. The pulse waveform after diffraction 
is compared with the incident waveform. The results converted from frequency domain by applying a inverse fast 
Fourier Transform (IFFT) will serve as sanity checks. 

Since waveform is distorted by the channel, we need to design transmit waveform in accordance with the different 
channels. Another contribution of the chapter is that we analyze die optimum transmit waveform based on derived 
time domain impulse response. 


3.2 Pulse Distortion and Its Impact on UWB System Design 

3.2J Pulse Distortion Caused by Wide Bandwidth 

Traditional communication theory are based on the assumption of distortionless propagation, which suggests that the 
receiving signal consist of a series of the replica of the transmitting signals. This assumption is valid when the signal 
bandwidth is not wide. However, for wideband signals, especially pulse based signals like UWB, this assumption 
might not hold anymore. As a simple example. Fig. 3.2 shows the waveform distortions for pulses with di ferent 
bandwidths after they pass through a 90-degree wedge (diffraction) channel. The incident pulse in Fig. 3.2(a) is 
a second order derivative of Gaussian waveform with a pulse width of 0.2ns, corresponding a 5.86 GHz I (MB 
bandwidth in the frequency domain, and the pulse in Fig. 3.2(a) is about 1 0ns wide, corresponding a 128 MHz 
10-dB bandwidth. For the convenience of waveform comparison, the amplitude of both incident waveform and the 
diffracted waveform have been normalized. It is evident that the shorter pulse (wider bandwidth) experienced more 
distortion after they passing through the same diffraction channel. When the signal bandwidth is not very wide (c.g„ 
the 128 MHz waveform in Fig. 3.2(a)), the distortion caused by channel is negligible. However, when the signal 
width is increasing, the distortion becomes noticeable. The details of the calculation of the diffraction waveform for 
a wedge diffraction will be given in the Section 3.5.1. 
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(a) Wide pulse diffraction (b) Short pulse diffraction 

Figure 3.2: Pulse distortion for pulses with different bandwidths. The incident pulse for (a) is about 10.2 ns wide, 
corresponding to a 10 dB bandwidth of 128 MHz, and for (b) is about 0.2 ns wide, corresponding to a 10 dB 
bandwidth of 5.86 GHz. Amplitudes were sealed to the same level to better compare the waveforms 

3.2.2 The Impact of Pulse Distortion on System Design 

Consider a communication system illustrated in Fig. 3.3. Here, p(t) is the transmitted pulse waveform, and h(I ) 
denotes the CIR. The received signal before q (t) can be written as 


r(t) — p{t) * h(t) + n(t ), (3.1) 

where n(t) is Additive White Gaussian Noise (AWGN) with a two-sided power spectral density of /Vi/2. Given 
s(i), a matched filter matched to s(-t) provides the maximum signal-to-noise power (SNR) ratio at its output. The 
matched filter here can be virtually decomposed into two filters, with the first filter matched to the CIR h(i) and 
the second matched to the pulse waveform p(t). As shown in Fig, 3.3, this decomposition can be done by setting 
d (£) — h(—t) and c 2 (f) = p(-/). conventionally channel is modeled as a tapped delay line. 

L 

k{t) = £>*(*■ r ') ° 2) 

[=1 

where L is the number of multipath components and 0 } and t ( arc its individual amplitudes and delays. The received 
signals are just a series of replicas of the transmitted signals, with different attenuations and time delays. 

In practice, the estimation of h(t) is generally a very difficult task, especially when there is per-pulse distortion 
for the wide band signals. The receiver suffers significant performance loss (as big as 4 dB), if the system is 
designed without the necessary compensation for pulse distortion [10], Along with system performance loss, pulse 
distortion can cause timing and synchronization error as well. It has been shown in (10] that the timing error caused 
by distortion due to cylinder diffraction can be much larger than the Cramer-Rao lower bound, thus is another 
fundamental issue that should be considered in system design. 
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Figure 3.3: General model for data transmission 


A new channel model that takes into account per-path distortion is then proposed in [15] 

L 

h(t) = £ hi(t) * 6{t - r,) (3.3) 

1 = i 

where distortion factor h{(t) includes both attenuation and distortion of the signal when pass through the channel. 

Note that waveform information can only be provided by time domain framework, which justifies our effort for the 
tune domain analysis. 


33 Time Domain Channel Modeling 


In this section we study several simple models that is mathematically tractable. For these models, the analytical 
results can often be obtained through exploring the physical mechanism behind it These analytical results provide 
us unique insight into the signal propagation in different practical environments. In this section, we are aimed to 
derive the time domain channel impulse response for these models. Wc will firstly start with its frequency domain 
model and then transform it to time domain via Laplace transform. 


3-3.1 Wedge Channel 


Rather than seeking the exact solution for Maxwell’s equations, Geometric Theory of Diffraction (GTD) uses gen¬ 
eralization of Fermat’s principle to find high-frequency harmonic Maxwell equation solutions for many complicated 
objects. The Uniform Theory of Diffraction (UTD) was proposed to complement GTD in the situations where GTD 
was invalid [4]. The original GTD/UTD was proposed in frequency domain and then was widely used in narrow- 
band radio propagation, e.g. s for path loss prediction. In 1990, time domain UTD was firstly derived by Vemttipong 
and Kouyoumjian in [5] by applying an inverse Laplace transform to the corresponding frequency domain results. 
TD UTD diffraction coefficient for a straight Perfect Electric Conduct (PEC) wedge can be expressed as [5] 




(3-4) 
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where 


and 


The value of X m in (3.6) is given by 
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(3-5) 


(3.6) 


(3.7) 


where /V 11 is the nearest integer solution of equation JbrniV* — 0 1 = ±7r, and L is the distance parameter determ ined 
by die incident wavefront. For plane wave incidence, we have I, — r sir r 2 /^, where r is the distance along the 
diffraction ray from the edge to the observation point. The rest of parameters used in (3.4)-(3.7) are described as 
follows: 


u s, h" soft and hard boundary conditions 

n the parameter describing wedge angle; wedge angle <fc\ = (2 — n)ff 

0o the angle between the edge and the incident ray 

<f> the angle between the wedge surface and the diffraction ray 

jp 

0 the angle between the wedge surface and the incident ray 
These parameters arc also illustrated in Fig, 3.1, which shows a special case when - tt/2. 


If the observation angle is not in the transition zone, we have the time domain diffraction coefficient for GTD (Keller) 
solution [5] 


- ——T ---tttt- (3.H) 

cos - — cast ——) cos £ — cos( ——) 

It should be noted that Veruttipong’s TD wedge diffraction coefficient given in (3,4) is a special case of the diffraction 
coefficient for curved wedge in [6], where an analytical signal representation for the transient fields has been used. 
The details for analytical TD-UTD will be given in the Section 3.5.1, 


D ah (t) = J -- 



sin 


n 


2 1 T17T Kill 0o 


3*3.2 Rectangular Building 

The first propagation model of interest is rectangular model. Rectangular is one of the typical shapes often encoun¬ 
tered in a practical propagation environment. Moreover, a lot of obstacles in the propagation environment can he 
approximated by a rectangular, with a limited degree of accuracy [7], ! his model has been widely studied in the past 
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Figure 3.4: UWB pulse diffracted by a rectangular building. Tx: Transmitter antenna, Rx: Receiver antenna. 
denotes the distance between Tx/Rx and the building. The size of the building is 2 a x b. The lower half plane plotted 
by dash line denotes the image of the building relative to the ground. 


decades. However, most of research is focused on frequency domain such as path loss prediction under different 
frequency band. In this chapter, we investigate its time domain characteristics such as pulse distortion introduced by 
diffraction in the edge of the building. Wc also analyze the optimum transmit waveform based on this time domain 
channel model. 

As illustrated in Fig. 3.4, we consider a simple UWB channel consists of a transmitter, a receiver, and a rectangular 
building that blocks the LOS transmission. In Fig. 3.4, dr denote the distance from the building to Tx, Rx, 
respectively. The size of the building is 2 a x b. The lower half plane plotted by dash line denotes the image of the 
building relative to the ground and is considered as perfectly absorbing. 

A good analysis of the path loss prediction in frequency domain can be found in [7] and the references therein. 
Generally, two approaches, Fresnel-Kirchhoff scalar approach and GTD approach considering a four-ray propagation 
system, are applied to study the effects of rectangular building diffraction. 


Propagation Model Based on Fresnef-Kirchhof Integral 


AIJ 1 Ycquency Domain Model 


Fresncl-Kirchhoff is one of the traditional approaches used to calculate the propagation loss due to the blocking of 
a rectangular building. We Define the channel transfer function as the ratio of the strength of the received signal to 
the transmitted signal. Mathematically transfer function can be expressed as H(jus) — Then the transfer 

function describes the rectangular diffraction can be expressed as [7]: 



(3.9) 


where A: is the wavenumber and d is the propagation distance. In our case wc have d — 4 + d T . In (3.9), variable A 
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and B arc the normalized half width and normalized height of the building, given by 


A = a* 


2 (dt + dr) 

A dtA 


B = b< 


2 (dt -T cl r ) 
A dtdr 


(3.10) 


where A is the wavelength. It should be noted that the validity of (3,9) is based on the condition Hi at the distances to 
and from the obstacle must be greater than its dimensions and the wavelength, he., 4, <4 >> u, b, A. 

A1 .2 Time Domain Model 


In this section, we will derive the time domain channel impulse response based on (3,9). We firstly work on one of 
the integrals and formulate them to a form ready for Laplace transform. 


it 2 , 


J exp(-J-y‘)r*/ 


r 

j ti 


exp i -( 


#> 2 } 


dy 



A* 


— : / exp(— z )dz 
kj Jo 


1 




where the complementary error function is defined as 


€T 


and 


2 r 00 

fc(x) = - 7 = exp {-t 2 )dt. , 

Jx 




y/& T a , 


Here s — jut , the propagation delay r a — a z {dt + d r )/(2d t d r c} % and c is the speed of light. 


(3.11) 


Similarly we have 


u 


7T 


/ exp{-j% 2 )f*/ 

Jo £ 


^(1 - j)(l -eTfc(y/W b )) 


where — tP(dt + d r )f(2d t d r c). 


Substitute (332) into (3.9) and after some manipulations wc can reformulate (3.9) as 

II{s) = 3 [1 _ (1 _ erfc(i/sf£)) (1 - erfc(^sr b ))\ e ST " 
where tq — d/c is the common time delay. 


(3.12) 


(3.13) 


Recall the Laplace transform pair 


er let Jin & _!=„(, - T .) 


7T ty/t-T a 


where the unit step function **(/) is defined as 


“ (f) - {i; 


t > o 
t < 0 


(3.14) 
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By applying inverse Laplace transform to (3-13), time domain impulse response can be obtained 


h{t) = <S(i-T 0 )* 


I / / u I r b \ f - r 4 ) it{t z T(,) ^ 

2 nt\\t-T a U T “ \t-T b b ) 2 n 2 - r Q ty/t - r b ) 


. (3.15) 


When a —► do, r a do, thus (3.15) reduces to 

h{t) = %S{t - r 0 ) * ^ ~ T b) - 0- 1 6) 

which is exactly the time domain formula for the well-known expression for knife edge diffraction [P.503, [17]]. 
This serves as our sanity check* 

Eq. (3.15) consists of three parts. The physical explanations for the three parts are straitforward. They represent the 
response from horizontal edge, vertical edge, and the mutual interaction between these two edges, respectively. 


Propagation Model Based on GTD 


Ray concept is involved when we attack the diffraction problem with GTD/UTD, The total received field strength is 
calculated by adding contributions from the various rays, which is also referred to as multipath components in the 
wireless communication. [7] provided a good GTD-based model in the frequency domain 


//CM 


+ ir' h (ju, <jf, <t>) Jkid ' 2 +d >>, 


(3.17) 


where the frequency domain GTD diffraction coefficient for half plane (zero-angle wedge) V^ h {jut f can be 


expressed as 

l 1 

cos ( n 2 R ) ^ COtS() 

The corresponding angle and the distance parameters in (3 AT) and (3,18) arc given by 


D sh (j^n\n) - 


4tv 


(3. IS) 


$ f = tan l (^) 

4>’ = 

rf] — df + (i 2 
d% — + IP 


i> = ip* + 7T + tan l {£) + tan 1 (£) 
4> = i// + 7T + ti mr* (g-) + tan~ l (£) 
4 = \fd% + a 2 
4 = sjdl + 6 2 


(3.19) 


The key term to transform (3.12) to time domain is the GTD diffraction coefficient, which here can be obtained by 
selling n — 2 and fa = n/2, from the GTD soltuion (3.8) for wedge diffraction in Section 3.3.1 


D sM {t, Si', il) 



(3.20) 
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Figure 3 J: UWB pulse diffracted by parallel plates. 

333 Two Parallel Half Plane Model 


Diffraction of a pulse or a transient wave (acoustic or electromagnetic) by a half plane is one of the early diffraction 
problems* Exact time domain solution of the CIR for such a channel has been obtained* In this paper, another 
case for the pulse diffracted by two parallel half planes will be studied. Again, we firstly bonow the existing exact 
frequency domain solution and then transform it to time domain by Laplace transform* It will be interesting to 
compare the field diffracted by two parallel half planes with that diffracted by a single half plane. 


In reality, a lot of environments can be approximately modeled as the problem of diffraction by two parallel planes. 
Fig. 3*5 illustrates an example. We assume the transmitting antenna is located far away from the planes and thus the 
incident wave could be viewed as a plane wave. Moreover, we assume the observation point is far from the planes 
so that the distance parameter p —> oo. We restrict that the incident angle is in the range of 0 < < n and the 

observation angle <p ^ ir ± which means 4> is not in the vicinity of the geometrical optics shadow boundary and 
the reflection boundary. The transfer function for such an channel would be [ 16] 


H(jv) = 


Xjr\ 2 sin 5^1 


cos co n&i 

cos [fca(siu <f> + sin 0 t} )] + 


/ i 1 _1_\ y' 

i/4nka \cos4> eotufri) (2m+l)^ 2 

c*'" / i | i \ V n 

VTrrfcfi \ cos ^ + cos <txy) ^ (2^ 

TTl— l 


cos |fca(sin <p — sin <£o)j - 
»s [fca{shi <f> + sin + O ( 


(3.21) 


where a is the distance between two planes. Note that another assumption made here is ku >> I, which can be 
satisfied cither if the distance a is large or the frequency is high. 


Taking a Laplace inversion (by setting s = -juf), the time domain impulse response can be obtained as 




5 (i — a/c (sin 4> + sin^u) — p/c) + 

6 (t + a/c (sin 4> + sin^j) - p/c) 
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where 


A 0 


\/2 1 sin^0sin|0o 


! sin k 




i i 



(3.23) 


Physical meaning for (3.22) is as follows. The first term in (3.22) consists of a superposition of edge waves from the 
two half planes, with the condition that each is excited by the incident pulse alone. The edge wave is characterized 
by a kernel oand can be explained with fraction calculus as “semi-integral”. 

/( 0 *^ = r ( i / 2 ) - £>“ 1/2 {/('■)} 

where l (x) is the Gamma function and the fractional integral of order 7 is defined by 


0- 7 {/(*)} - p^y fj - cr 'nm 

The second and the third terms are infinite sums, which represent the successive mutual interactions between the 
two half planes. These two terms are characterized by a kernel of ti(£) and we have 

m * n(t) - f 

The coefficient OT (iin)^ decreases rapidly as the value of m is increased. Practically, computation of 

infinite sums in (3.22) requires but few terms to obtain reasonable accuracy. This will be verified in our numerical 
results in the Section 3.5. 

It should be noted that if both the transmitter antenna and the receiver antenna are located in the cannon between 
two half planes, then there is no diffraction phenomena and multiple reflections dominate the received signal. Time 
domain C 1 R for such a channel has been studied in [18]. 


3.4 Time Domain Optimum Transmit Waveform Design 

Section 3.2.2 shows the optimum receiver should be a matched filter matched to the waveform distorted by the 
channel, not the transmitted waveform itself. However, from system point of view, such a matched filter alone is 
not enough. We can further maximize SNR at the receiver side by carefully designing the transmitted waveform 
[ 6 ]. The primary difference between the problem in this paper with the standard matched filter problem is the pulse 
distortion caused by channel. We will show later that if channel is distortionless, then the SNR at the receiver side is 
independent of transmitted waveform and then there is no optimum transmitted waveform. 


3.4J Optimum Transmit Waveform for Distorted Channel 

The transmitted waveform optimization problem can be stated as follows. Given the channel impulse response h{t) 
and fixed transmitted power F\ (average power), we wish to achieve the maximum SNR at the receiver side by joint 
design the transmitted waveform and a good receiver. This problem has been discussed in [6 [ for radar detection. 
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Assuming the transmuted pulse p(f) (to be optimized) is confined to the symmetric time interval \—Tf 2, T}2). The 
energy of transmitted pulse is then 



(3.24) 


It follows from detection theory that the best receiver is still a matched filter matched to the signal component in 
r(£), the SNR after such a matched filter is given as. 


SNR = 2E w /Nq* (3*25) 

where E y | 7 ^ 2 |p(t) * h(t}\ 2 dt is the received signal energy. The problem is then reduced to find the optimum 
p(i) such that Ey is maximized, under the constraint of fixed Ep. 

It has been shown in [13] (p, 125) and [6] that the optimum p(t) can be obtained by solving the following homoge¬ 
neous Fredholm integral equation 


fTf 2 

Pn <Pn{t) - I «(* “ T)4>n(r) dr , (3.26) 

J-T/2 

and lei p(t) - where ) is the eigenfunctions corresponding to the maximum eigenvalue /n and the time 

reversal operator /c(i) is defined as *c(f) = h(t) * h( —t). Witliout loss of generality, it has been assumed that 
po > p| > p2* * - * ■ Pulse waveform <fo(f) reproduces itself, sealed by a constant p,, when convolved with the time 
reversal operator over the interval [-Tf2 r T /2], 

With optimum p{() and c(f), wc achieve maximum SNR 

SNR = 2tioE p /N 0 . (3.27) 

It is worth noting that to meet the requirement of optimum information transmission might result in another transmit 
waveform. Maximization SNR requires the transmitter to put as much as possible energy to the frequency with high 
transmission capability through the channel, under the fixed transmitted energy constraints, while ignore the other 
frequency modes with small transmission capability. It is possible that these smaller modes contain a significant 
quantity of information useful for the system. We consider an extreme case here. When T —* oc, it has been shown 
in [15] (p. 360) that all the energy should be placed into one single frequency fnax in the operation band for which 
//(/) is the greatest. This leads to a transmitted waveform of a sine wave at jmax* However, it should be noted 
that a signal of arbitrarily narrow bandwidth conveys information at a zero rate. 1 herefore, trade-off should be made 
between SNR and information transmission. 

To illustrate the principle, we study the optimum transmit waveform for the following channels. Consider practically 
all the communication systems are band limited, we artificially impose a band pass filter in the system. The transfer 
function of this filter Hp(f) is defined as 



/i < I/I < h 

otherwise 


where /( and f 2 arc the low and high frequency of the filler, respectively. Its time domain impulse response can be 
expressed as 

hp(t) ~ 2f2Sinc(2f2t) — 2 fisinc( 2 fit) 
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The autocorrelation of the impulse response of this filter can be written as 

hp{t ) * — [2/251710(2/2^) — 2 /j $inc{ 2 f\t)\ * (2/2s£tic(”2/2£) — 2 f\si 7 ic(— 2 /ji)] 

= 2 /25inc{2 / 2 O + 2/] sinc(2f\ t ) — 4 f\ sinc(2f\ t) 

= h F (t) (3.28) 

Apparently autocorrelation of the function hf'(t) reproduces itself. 

ft should be noted that for cognitive radio purpose, the frequency response of this filter could be in an arbitrary 
passband (e.g. t kih element subband) located under the FCC spectral mask (3.1 GHz -10.6 GHz ) [19]. 

Assuming the CIR of a specific channel then the overall of the CIR will be h(t) ~ hp(t) * h c (f). 


Wedge Channel 


Consider a case (hat light of sight (LOS) is blocked by a PEC wedge, as illustrated in Fig. 3.1, the CIR of such a 
channel without consideration of bandwidth can be expressed as 


ml 


where 


/l c (f) - ri7r sin/ift 

= <7(n, ft, 


* fc sin^ 

C(n,f3o f 4>i$ “ J 9 ’ -dfe 
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COS 


£ — cos( ) cos ^ — cos( - ) 

n ' n f n ^ n > J 


rr T 




is a constant whose value is depending on the environment. 

Note that u{t)j\ft is an even function. The autocorrelation of the above CIR is 

. ,, w , , rt _ n -i "(0 . u (0 

M) M )_ Vt yft 

Then the autocorrelation of the overall CIR for a wedge channel is written as 

k (£) — h c {t) * h c (-t) * hp(t) * h F (~t) 

- c 2 -r 2 (i/2)-ir l {h F (t)*h F {-t)} 


- C 2 tt ■ f [ 2 / 2 sine(2/^f) - 2}\sinc{2}\i)\ 

J—OO 


dt 


= C* • [Si (2 jt/i t) - S*(2jt / 2 i)] + CV/2 


(3.29) 


(3.30) 


where sine integral (unction Si(x) is defined as 


Si{ 


« - II t 5 * 


Note that Si(x) can be implemented with a Matlab built-in function, substitute (3.30) into (3.26) and solve (3.26) 
for </>(/) numerically wc can obtain the optimum waveform for wedge channels. 
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Rectangular Channel 

The time domain CIR for such a channel has been derived in Section 3.3.2, It follows from Eq, (3.15) that the 
CIR consists of three parts. Practically, the contribution from the third part (mutual interaction between two edge 
waves) in (3.15) is very weak and could be ignored. We will also verify the validation of this assumption through 
our numerical results subsequently. Then the CIR can be approximated by the first two parts 



{t + T a ) \ft 


U(*) * <S(f- - T n ) + 


1 

(t + Tfj) v7 


u(t) * 6(t - T h ) 


(331) 


Note in (331) we have ignored the common time delay tq, for simplicity. 

For the term —-— 1 4 most of the energy is focused at the singularity where t — 0, thus wc have 7 -- — r 

(i+Vklv't’ a j (t+r„, t )vu 

— Then (331) can be further approximated as 


h(t ) * 6(t - r a ) + Cb^~ * S(t - Tb) 


ft 


y/t 


( 332 ) 


where = 237^ and Cb = 53^' 

The time reversal operator n(f) for such a rectangular channel can be expressed as 

<t) * (C 2 + Cl) ■ [Si( 27 r/, t) - Sifaht)] 


+C a C b ■ [S*(2jr/,t) - Si{2wf 2 t)] * \6{t - t b + r„) + 6{t + r n - r h )] + {C a + C fc ) 2 ir 2 /2 


( 333 ) 


Two Half Planes Waveguide Channel 


As slated tn Section 3.3.3, the half planes waveguide model is a little more complicated than the other two models 
studied above, due to the exist of multiple replicas (caused by waveguide phenomena). However, wc can simplify 
the problem by focus on the main lobe of the n(t) instead of the whole *(£}. The main lobe of k(1) for two half 
planes channel can be approximated as 


K(t) » 2Ag ■ [Si(27r/it) - Si(2nf 2 t)\ 


+24 ESU m- 3 ) - (Si(2*ht) - Si(2nf 2 t)) + & (^M _ “=1*6*1)] 


(3.34) 


For the case of high frequency where f\ —* 00 and / 2 —* 00 , interestingly, k(1) is reduced to a form of 


k(0 SB C f • [(Si&wht) - Si( 2tt/ 2 *))j 

where C f = £m-i is a constant depending on the environment parameters such as the location 

of the receive and transmit anlennas(^, 0 T ^), and the distance between the two half planes (a). Therefore, under 
these assumptions, the optimum transmit waveform for the two half planes waveguide channel will be the same as a 
wedge channel, which has been solved above. 
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3.4.2 Numerically Computing Eigenfunctions 

Except for some classical case such as band limited signals [13, 6], the eigenfunctions and eigenvalues are gener¬ 
ally difficult to be solved in a dosed form. In this section we are interested in computing the eigenfunctions and 
eigenvalues numerically. The idea is to discretize the function and borrow the powerful matrix theory to solve the 
problem. We firstly discretize the autocorrelation function ) with sufficient low sampling rate such as to include 
all the significant channel frequencies. 

As stated in the previous section, the optimal pulse waveform (maximizing SNR) is a solution to the integral function 

fT/2 

Pn4>n(t) = / - T)<t>n(r) & ( 335 ) 

J T/2 

Note that the TR operator K(t) is an even function,, i-e., K(t - r) — k{t -f), the above equation can be reformulated 
as 

f T/2 

= / *(T - t)4>n(r) dr (336) 

J -T/2 


r-T/2 

fLn * <£n(0 = / K(t - T)0 n (l) dr n=l,...,oo, (337) 

JT/2 
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M~N + 1) 

* 

*■ 
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#■ » 

■ 

k{0) 
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M-H) 

M-M + i) 

♦ 

- At 

(3.38) 


<M*) 


i# 

k(0) 

* 

k(2 JV) 


0 n(/V) 




We then transform the integral function to a matrix problem. To find the solution of the integral function (335) is 
then equivalent to finding the eigen function of matrix n(t)(2N -h i x 2N -f- I) in (3.38). Several remark notes can 
be drawn from (338) 

• The optimum waveform obtained via the numerical solution introduced above is not necessarily the unique 
one; though the maximum eigenvalue is unique, the corresponding eigenfunction might not. 

• If there is no pulse distortion (e.g., free space transmission), the time reversal operator *(f) — q $6(t - To), 
where qq and tq are the constants depending on the distance between the transmitter and the receiver. For this 
special case, there is no optimum transmitted waveform since any waveform will reproduce itself after being 
convoluted with the time reversal operator 

• We can improve computation accuracy by decreasing sampling interval At. However, the computation com¬ 
plexity is increased correspondingly* 
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It is also evident that the availability of time domain OR h(t) (and thus its time reversal operator «(/)) is critical 
to computing the optimum waveform The following will show how to get the time domain OR for several 

canonical environments. 


3.5 Numerical Results 

3.5.1 Calculation of Singularity in the Impulse Response of Wedge Diffraction 

A singularity appears in the impulse response for most of channels where diffraction phenomena is involved, e.g., 
in Eq. (3.6) (wedge diffraction), when t —*> 0, then IP> h (t) —> oo, thus in Eq. (3.4), there is a singularity at 
t = 0. How to numerically calculate this singularity is not a trivial issue. Since the singularity is exactly where 
the energy concentrates and also defines the type of the pulse distortion, it can not be directly removed using the 
time windowing. From an Electromagnetic point of view, this singularity at the diffracted wavefront is essential. 
An integral'differential (Int-Diff) approach has been proposed and proved to be an efficient way to remove this 
discontinuity at t = 0 during the calculation [3,10], In this section, we will further verify this approach by comparing 
it with the well known ATT approach [8], 


Int-Diff Approach 


Mathematically, the principle of this Int-Diff approach can be described as follows. 



(3.39) 


hence, we have 


P(0 * 9(t) = 



(3.40) 


where g(t) has a singularity in the time domain, The integration J,j g{t)dt removes the singularity and can often be 
obtained in a dosed form. The reason why we call this approach Int-Diff approach is because additional integration 
and differentiation operation arc involved to get the convolution result. 


Int-Diff approach provides a very efficient way to do the numerical convolution for the case when there is a singularity 
in the functions to be convolved. This approach is specially important when we deal with the diffraction problem in 
the lime domain because there is always a singularity in the diffraction coefficient due to the nature of diffraction. 

Again, we consider wedge diffraction case as an example. Notice that the singularity in (3.4) is solely caused by term 
F(X m J). We will focus on this term first. The rest of work is just linear summation. According to int-diff approach, 
firstly we will need to integrate Happily, we have a closed form for the integration of f), which 

can be expressed as 


(3.41) 
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Let D*' h (t) denotes the integration of the wedge diffraction coefficient I7 h (t), then wc have 


D tJl (t) = 


-i 


2fiv^27T sin(3o 


Y K^ h F{X m ,t) 


(3.42) 


Note that singularity in has been removed in the new integration expression D* th (t). The biter can be used 

to do the convolution, followed by a differentiation operation as a compensation. 


Analytic Time Transform (ATT) 

AT I is another classic technique used to perform an efficient convolution with a broad class of excitation pulse 
functions ini-diff. In this paper, this classic approach will be briefly introduced and serve as a sanity cheek for the 
int-diff approach. 

The second version ATT in [8] is defined as 


/w 



far Im(t) > 0 


(3.43) 


The "+** sign over / means an analytical signal in a complex-time domain. Let p(t\ and h(t) denote the incident 
pulse and the impulse response of the system, respectively, then we have. 


p( 0 * Mf) = \ 


P (t)* h (f) 


(3.44) 


+ 


Notice thal the singularity on real time axis (Im(t) — 0) in ft (f) is only a branch point singularity and not a pole. 

Therefore, the difficulty in calculation the singularity in ft(£) has been removed by ATI in ft {/), Moreover, if wc 
choose the incident pulse such that in the frequency domain p(cj) can he modeled using exponential function as 


N 


p{ui) = A n e anUt for oj > 0 


(3-45) 


n- f 


or equivalently, in time domain 


N 


P (f) = X] 71 " 6 ( t + j°‘n) 


(3.46) 


n —1 


then the convolution can be evaluated in closed form as 


p(£) * h{t) — Rc 


N + 

^ ^ A n ft (t H- jot n ) 


.n — I 


(3,47) 


where the shifting properly of the analytic delta function 


l - 6 {t+ja n )*h{t) =4 (t+jan) 


(3.48) 


has been used. 
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Figure 3.6: Comparison between the approachs of Int-DifT and ATT, In (a) amplitudes have been normalized to the 
incident pulse to compare the power of the incident pulse and the received pulse. In (b), the incident pulse has been 
scaled to the same amplitude of the received pulse to compare their waveforms 


To verify our result with int-diff approach, we use the same incident pulse as used in [8], Mathematically, the 
frequency domain of the incident pulse p(w) is defined by 

p{uj) = C 0 { 1 - e~* T ) p 'e‘“ F * T (3.49) 


where I' — j^FT^ 71 

ent requirements. The peak of the p(u>) is normal i/xd by choosing ( i — ( 1 * Time domain 
waveform of the incident pulse can be conveniently obtained by 


( l [ p y hi ), }<■ represents the center frequency, which can be adjusted easily according to dilTcr- 


tt 
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(3.50) 


where = 

Let I) (0 denote the analytic impulse response of wedge diffraction, then using the convolution property in Table 
2.1 in [8], the receiving pulse will be 



Re 

Re 


{ p (t)* b (0 

J 

b {I + j(n + R 2 )T) 


(3.51) 


Instead of integration, ATT approach remove the singularity by performing an analytic transform. The analytic 
transform of an time domain function is not so straightforward, especially when the function itseh is complicated. 
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According to [8], the analytic wedge diffraction coefficient I) (/) can be expressed as 
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-1 


2n\/27rsmdt|j 


) ' m F (^mi 0 


m—1 


(3.52) 


where 
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+ yjjx m /c) 


(3.53) 


and the rest of the parameters are the same with that of TD-UTD, presented in section II. 

Fig. 3.6 shows the received waveforms, calculated with ATT and Int-Diff approach, respectively. The incident 
waveform, marked as red slash curve, is also plotted, as a reference. Firstly, it can be easily observed that the 
received waveforms, calculated with different approaches, are in excellent agreement. Secondly, it can also be 
observed from Fig. 3.6(a) that the diffracted pulse is much weaker than the incident pulse, due to the huge energy 
loss caused by diffraction. To better compare the waveform between incident pulse and the diffracted pulse, we 
manually scaled their amplitudes to the same value and their waveforms are shown in Fig. 3.6(b). It is evident 
that the diffracted pulse waveform is much different from the incident one, indicating a pulse distortion caused by 
diffraction. Finally, the amplitude information (signal strength) can be predicted by traditional frequency domain 
UTD, the pulse waveform information that presented in this paper, however, can only be provided by lime domain 
formulations. 


The parameters used in our simulation are as follows: For the incident pulse, H = 2, P 2 — 1, and / c = 2 GHz. The 
parameters for the wedge are: $1 — rr/2, n — 3/2, 4> = 57t/6, <t> ~ jt/ 3, and the so ^ boundary condition has been 
considered. 


3,5.2 Rectangular Building 

For all the figures (unless otherwise specified) in this paper, the waveform labeled “TD" is obtained by convolution 
of the incident pulse p(() with the time domain CIRs derived in this paper, and the waveform labeled “FD Hi IT’ is 
obtained by applying inverse Fast Fourier Transform (FFT) to the product of P(f)H(f) 7 where P(f) and //(/) are 
Fourier Transform of p(f) and channel transfer function, respectively. We have verified our time domain derivation 
by comparing the "TD” result and “FD+IFFT” result. 

In the following numerical results, we choose the second order derivative of Gaussian pulse as the incident pulse, 
which mathematically is defined as 



(3.54) 


where w is the parameter controlling the width of the pulse (and therefore the frequency bandwidth of the transmit 
signal) and t c is time shifi to pul the pulse in the middle of the window. 

Fig. 3.7 shows a comparison between GTD and Fresnel-Kirchhof models. Here, we set w - t). I ns and — 0.3 
ns. The environment parameters used in the simulation are as follows: 4 = 1000m, dr — 900m, a 15m, and 
b — 10m. 

Several observations can be made from Fig. 3.7. 
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(c) Comparison Kiichhoff ami GTD results 

Figure 3*7: Simulated pulse waveforms after a second order derivative oi Gaussian pulse pass through a rectangular 
channel described in Fig, 3,4 


• For both GTD model and Fresnel-Kirchhoi Model, time domain results show very good agreement with the 
results transformed from the frequency domain via IFFT (labeled “FD+IFFT”). This verifies our lime domain 
derivation. 

• The diffracted waveforms based on GTD and Frcsnel-Kirchhof agree very well to each othci\ except that 
the Fresnel-Kirchhof result has an additional small waveform response. This extra “baby response’' can be 
identified from its time domain formula to be the mutual interaction (coupling) term between two edges, which 
has not been taken into account by the GTD-based model. 

Again, for the convenience of comparison, common propagation delays has been neglected. 


3.53 Two Parallel Half Planes Model 

As mentioned in Section 3.3.3, when a UWB pulse is incident on the parallel half plane channel, the received signal 
consist of three parts, with each part including two same (but separated in the time) signal components. The relative 
signal strengths of these three parts depend on the value of the coefficient and Aj, which in turn depend on the 
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Figure 3,8: UWB pulse diffracted by a parallel plane waveguide 


environment and the location of the antenna. Fig, 3,8 shows the received signal when a second order derivative of 
Gaussian pulse is incident on the parallel half plane channel. The incident pulse width parameter w — 2 ns. The 
geometrical parameters for parallel planes are as follows: p — 20m, <j> = 8/15tt, — 5/6rr, a — 2m. 

It can be seen from i ig. 3.8 that the waveforms of the received pulses are not second order derivative Gaussian 
anymore. It follows from Section 3.3.3 that they are actually “semi-integral” forms of the second order derivative 
Gaussian pulse. Again, the agreement between the waveform obtained by time domain convolution (labeled “TD”) 
and frequency domain result (labeled “FD+IFFT”) shows our time domain derivation is correct. 

Although Eq. (3.22) shows infinite terms, we can see from Fig. 3.8 that these sums converge very fast. 


3,5,4 Optimum Transmit Waveform Design 


Fig. 3.9 shows the optimum transmit waveforms for different channels. These waveforms have been obtained 
by substituting /c(£) of corresponding channels into (3.38) and numerically computing <p{t) associated with the 
maximum eigenvalue. For the convenience of waveform comparison, the energy of the computed optimum transmit 
waveforms for different channels have been normalized. The parameters used in the simulation are as follows: 
T — 2ns y ft = 3 GUz y = 5 GHz. Moreover, in each model, the environmental parameters we have used here 
are the same as those we have used for the computation of distortion waveforms in this section. It has been show n 
in Section 3.4.1 that the optimum transmit waveform for the two parallel half planes channel is approximately the 
same as wedge channel, under some circumstances. 

For the purpose of maximizing the SNR at the receiver side, all the energy should be sent through the optimum 
transmit waveform, which is the eigenfunction associated with the maximum eigenvalue of #c(t). For the purpose 
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Figure 3.9: Optimum transmit waveform design for bandlimited UWB signal passing through different channels 

of flexibility, we are interested in transmission of a linear combination of the eigenfunctions whose eigenvalues 
are close to the maximum one. Fig. 3.10 shows the first six eigen waveforms for the passband channel. Their 
corresponding eigenvalues are 1,0000, 0.9999, 0.9979, 0.9976, 0.9621, and 0.9599. We can also transmit these 
waveforms parallel to increase the data rate. Note that the eigenfunctions for the low pass channel {(0, jfj) is knows 
as prolate spheroidal wave functions [20], and has been proposed as a novel UWB pulse shaping technique in [21 ]. 


3.6 Conclusion 


For wideband signals such as UWB, by principle, time domain frame work appears to be a better frame work to study 
its propagation. The time domain impulse responses for a series of scattering environments have been derived for 
the first time. Specifically, several electromagnetic communication fundamental channels, including perfect conduct 
electromagnetic wedge, large rectangular plates, and parallel planes waveguide, arc studied. The receiving signal, 
calculated through convolution in the lime domain impulse response, can provide waveform information, as well 
as the amplitude information (Path loss). The pulse waveforms passing through these channels arc compared with 
their original waveforms. The comparison show that for wideband signals, a transmit pulse could be extremely 
distorted by diffraction mechanism. The results converted from frequency domain by applying a inverse fast Fourier 
Transform (IFFT) serves as sanity checks for our time domain derivation. Given the lime domain channel impulse 
response, we also study the optimum transmit waveform such that the signal to noise ratio (SNR) at the receiver side 
is maximized, under the constraint of fixed transmitted power. 
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Figure 3JO; The first six eigen waveforms for a passband channel 
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Chapter 4 


Optimum Transmission Waveform in RF 
Harsh Environments 

4.1 Introduction 


Communications and networking in confined metal environments are critical to applications such as intra-ship, intra- 
vehicle, intra-engine, manufacturing plants, assembly lines, nuclear plants, body area network sensors surrounded 
by vehicles and tanks, etc. Narrow band wireless technologies have proved ineffective in these environments, due 
to resonance caused by the metal walls. The physical mechanisms of short pulse propagation in a confined metal 
environment are studied by Felsen[l ], from a transient radar cross section. An Ultra-wideband (UWB) system [2]- 
[5] is proposed here to resolve the resonance into many time-resolvable pulses—corresponding to extremely rich 
multipath. The energy of these pulses is difficult to collect using traditional RAKE or multicanier technologies. As 
an alternative, a powerful time reversal channel-matching technique is suggested for this purpose. 

Periodically modulated radar returns from jet engines can be modeled as a chirp [6]. This view provides some insight 
into our problem, with measured channel pulse response from a pulse train as long as 800ns. 

Chirp UWB systems [7]—accepted into IEEE 802.15.4a in Dec. 2Q0fV have some favorable advantages such as 
low cost, low transmission power (long range), low battery power, anti-jamming, easy synchronization, robustness 
to Doppler frequency shift (mobility), etc. Analog SAW filters can be employed to perform receiver matched filter. 

This paper is motivated by die UWB applications in the harsh environment, by adapting the existing standards. 
Time reversal[9][I0]—see [8] for a review—has been used to compress the channel response into a sharp main 
lobe. Similarly, chirp matched filter results in pulse compression. We design the time reversal chirp system from 
the classical modal modulation and give some relative performances in different situations. Finally, the spectral 
efficiencies of different transmission schemes are compared to benchmark the time reversal chirp scheme. 
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Figure 4.1: Rectangular metal cavity used for propagation experiments. 

4.2 Pulse Propagation in Rectangular Metal Cavity 


The measurement is performed using frequency domain technique to analyze the characteristics of the UWB channel 
in the confined metal environment. The rectangular confined metal cavity in ig. 4.1 measures 16 feet by 8 feet by 
8 feet. The materials of the walls are aluminum. 


The sounding is carried out by sweeping a set of narrowband sinusoid signals (tones) through a wide frequency 
band. The data of channel frequency response is collected using a Vector Network Analyzer (VNA)—Agilent 
N5230A (300kHz - 13.5GIIz). The VNA sweeps from 3GHz to 10GHz using 7001 points with frequency step 
of IMIIz . The power of each tone is IQdBm, and the same for all the measurements reported in this paper. The 
number of averaging is 128, which is used to enhance the accuracy of the measurement. Meanwhile, the cables and 
the connectors are calibrated before measurement to compensate for the frequency dependent losses. 

The measurement is done for line of sight (LOS) situation. The distance between transmitter antenna and receiver 
antenna varies from 0.5m to 4m with distance step of 0.5m and both antennas are placed at 1.35m above the ground. 
The transmitter antenna is fixed, and the receiver antenna is moved to different locations along the middle line of 
one side of the rectangular metal cavity. 

After the channel response is measured in the frequency domain (Fig. 4.2), zero padding is done between 01 Iz and 
3GHz, conjugate reflection between -lOGIIz and GHz and inverse fast Fourier transform (IFFT) are used to obtain 
the channel pulse response h(t) in Fig. 43 and Fig. 23(d). Channel pulse response in the office environment 
(Fig, 23(d)) is very different from that in the rectangular metal cavity (Fig. 43), for the transceiver distance of 4m, 
Channel pulse response in the metal cavity is tong and strong for the transceiver distance of 4m (Fig, 43). 


The energy of the channel is 



(43) 


where Th is the time duration of the channel. The autocorrelation of the channel is shown in i : ig. 4.4, defined as 


Rtoit) = h(t) * h(T h - i). 


(4.2) 
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Figure 4.2: Measuring frequency response //(/) in a rectangular metal cavity. A total of 7001 equal frequency tones 
arc sent from 3GHz to 10GHz (with a step of I MHz). 
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Figure 4.3: Channel pulse propagation h{t) in a rectangular metal cavity. The pulse has 7001 equal frequency 
components (with a step of 1MHz) from 3GHz to 10GHz. 


It is interesting to notice that the total energy of the channel pulse response in a confined metal cavity drops very 
slowly as compared to office and hallway environments, as shown in Fig. 4 + 5 + 


4.3 Time Reversal Transmission with Chirp Signaling 

43.1 Modal Modulation 

The starting point of our system formulation is from the well known result of Holsinger (1964) [13], regarding 
optimum detectors for a linear time-invariant (LTI) channel with memory, as shown in Fig. 4.6.1: x(t) is the input to 
the LTI channel with channel pulse response h(t) mentioned before. The symbol duration is 1$, and the waveform 
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Figure 4.4: Autocorrelation of channel pulse response /) — h(t) * h( i' h — t)/E h in a rectangular metal cavity. 
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Figure 4.5: The total energy of the channel pulse response in a confined metal cavity decays very slowly, as compared 
with office and hallway. 


duration of x(t) is T. The output waveform is 


r 



x (r)h (t — r) dr , 0 < t < T\ 


(4.3) 


for the observation Lime [0 t T|]. White Gaussian noise is assumed. According to the channel pulse response, the 
observation time and waveform duration, a symmetric function K(t, s) can be defined as 


7*"i 

K(t,s)= f h(cr — t)h(a — s)da, 0 < t,s <1 (4.4) 

Jo 

Furthermore, we can define a possibly infinite-size set of orthonormal eigenfunctions, ^(/), that satisfy the relation 
[13] 

An • <Pn(t) = / K(t, T)<p n {r)dr, 0 <t<T 

Jo 


(4.5) 
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Figure 4.6: Modal modulation is optimum for the LTI channel with memory* When a self-reproducing eigenfunction 
ipn(t) is transmitted, after the channel matched filter h(min(Th/fi) - t) y the same waveform <f n (t) (within a gain 
factor A ri ) is obtained. The optimal detection waveform is matched to the <#*(*). 



I >TtT t 


Figure 4.7: Time reversal modal modulation is proposed to simplify the transceiver complexity. The channel matched 
filter h(T h — t) is moved to the transmitter side. The self-reproducing <th(t) can be defined, similarly to the modal 
modulation illustrated in Fig, 4.6* 


A special form of (4,5) is given by [14]. That is a function that when convolved with K (/, s) over the interval 
[0,7'] reproduces itself, scaled by a constant called eigenvalue* 

£ n {t ) n — 1, 2,.,., oo can be seen as a set of basis functions, so x{f) can be written as 

oo 

X (f) = ^ X n if n (t), 0 < t < T (4*6) 

n =0 


and 


x n — f x (t) <pn n — 1 T 2,.. 

Jo 

so the energy of the transmitted signal is 


., oo 


& 


= f i*(*)i : 

Jo 


dt 


CO 


= S*n 


n=l 


H.7) 


(4.8) 

(4.9) 


From detection point of view, after matched filter, the output signal to noise ratio (SNR) is 2 Ej IS}, where 

<x 

E = Y1 X n xl 

n= 1 


(4.10) 
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and /V ( |/2 is the double sided noise spectra] density. So the probability ol error for the detector is given by 





(4.11) 


In order to minimize P e , SNR should be maximized. Since by conversion A| > A% > A 3 > * • •, if the energy of 
x{i) is fixed, maximizing E means x(t) - xi^i(f) which is the optimal waveform for detection. So, the maximum 
SNK can be expressed as. 


SNR 


max 


2 


2Xixi 

“NT 

^max i‘ J t 

Nn 


(4.12) 

(4.13) 


From capacity point of view, the modal modulation system have many parallel orthogonal channels and each channel 
corresponds to one eigenfunction Wc can use the classical water filling scheme to calculation coefficients 

x„ n — L, 2 ,..., oo based on the energy constraint of transmitted signal x(/). 

The observation time 1\ and the symbol duration Tg are very important for the system design. If 7] < T + 7^ T then 
l'$ should be equal to or greater than max (7j + min(Tj,7 a)/T + 7k). And if T\ > 7 + 7 k, then K (f, s) reduces 
to K (t — s) and Vg should be equal to or greater than 7 H- 27k. In these situations, the generalized Nyquist criterion 
can be satisfied, so wc do not need to consider IS! when using modal modulation. 


4.3,2 Time Reversal Modal Modulation 


Time reversal modal modulation is proposed to simplify the transceiver complexity. The h(% - t ), lime reversal 
filter is moved to the transmitter side. Similarly to the modal modulation illustrated in Fig, 4,6, the self-reproducing 
eigenfunctions 4^(1) with eigenvalues ^ illustrated in Fig. 4.7 can be defined as 


where 


K n • 4>n{t) = f Hhh(t - r)4> n (T)dT, 0<t< I 

J o 

Rhh(t) = fi{t) 


(4.14) 

(4.15) 


Notice that when T\ > T + 7k and Tg > T + 27k, moving the channel matched filter ft(7k - t) in Fig, 4.6 to the 
transmitter side will not affect the optimality of the system. In other words, the systems in Fig, 4.6 and Fig. 4.7 
arc equivalent in the above situation. In Appendix B, we give the proof that when the optimal waveform is used for 
detection, the maximum SNIt after matched filter is equal to 2^ nax f^/jVo similarly as shown in Eq. 4.13. 


4J3 Time Reversal Chirp System 

Time reversal modal modulation gives us a lot of insights to design UWB system. The main work shifts from 
how to handle IS1 to how to design the transmitted waveform. We can calculate the orUionormal eigenfunctions 
4>n(t) n “ I> 2,.,,, oo according to the autocorrelation of the memory channel. Each eigenfunction can be seen as 
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Figure 4.8; Output waveform of the chirp matched filter without time reversal, defined as /(*.(/) in (4.17). 



Figure 4.9: Output waveform of the chirp matched filter with time reversal. 

a waveform and any linear combinations of these waveforms can be used in different situations and requirements as 
well. 

It is well known that chirp signal has the sharp autocorrelation characteristic and it is widely used in the radar field. 
We can easily use chirp signal as the waveform instead of the eigenfunction to get the time reversal chirp system. 

Chirp signal is defined as 

= 00 .( 2 . 

where /«, 1\ fi are center frequency, signal duration and sweep rale, respectively. The amplitude spectrum of a chirp 
signal approaches to a band pass filler as the time bandwidth product of it becomes large. 

The frequency of a chirp sweeps continuously, and its duration can be changed independently. ! lie bandwidth B of 
a chirp signal is B — ^ * 1\ Signal with +/i is an up chirp while —fi is a down chirp. The pair is matched to each 



(fut ± ^ j ),0<t<T 


(4.16) 

















58 


CHAPTER 4. OPTIMUM TRANSMISSION WAVEFORM IN RF HARSH ENVIRONMENTS 



Figure 4.10: Simulation result for time reversal chirp system in the rectangular metal cavity. 


other. The output of the matched filtering is Rcc( f ) = Oi(0 * Ch('T ~ t )—implemented using surface acoustic wave 
(SAW) devices—expressed by [19] 


= v^cos (2tv f 0 (t - '/’)) 


sin (t if.it (7' — \t — 'f\)} 
2 nfit 


0 <t< 27 


(4.17) 


which is illustrated by Fig. 5.3, 


if time reversal is not considered, after the chirp matched filter, a compressed pulse that has a short duration of 2 / B 
is generated. The peak in Fig. 5.3 is \/TB. As a result of matched filter, a fiat envelope pulse with duration T is 
compressed into a much shorter waveform of duration 2//?, If time reversal is used, the output waveform of the 
chirp matched filter is shown in Fig. 4.9. 


4.5A System Performance Comparison 

The system performance comparison between time reversal modal modulation system and time reversal chirp system 
will be considered here when 7s > T 4- 27V 


From Detection Point of View 

From detection point of view, SNR after matched filler in the receiver side is the metric to determine the system 
performance. 

The maximum SNR for Lime reversal modal modulation system can be shown in Eq. 4.108. 

For lime reversal chirp system, if chirp signal o, (f) is used, then, 

00 

ch w - Cn4>n w. °^ t ^ T 

n—l 


(4.18) 
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where 


f-TL 


-f 


c h (i) <f> n (i t)dt , n = 1^2, — ,oo 


and <f>n (t) is from Eq. 4.14. 

The energy of the transmitted signal is 
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In the receiver side, after matched filter, SNR is 
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(4.19) 
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So, ihc maximum SNR gap between time reversal modal modulation system and lime reversal chirp system is 


9SNR 


OO 

n-0 

Ef 

oo 2 
^max X^ c n 


oo 


E 

n—i 


(4.31) 


(4.32) 


If 7' goes to infinity, Cn can be treated equally and Eq, 4.32 is reduced to 


SSISIl = 


(max ([/-/ (/)!})' 


which can be treated as the peak to average ratio from perspective of frequency domain. 


(433) 


When T is equal to Ins, 10 ns, 100ns and goes to infinity, the maximum SNR gap will become 2.5dB, 4.5dB, 8.4dB 
and !2.7dB, respectively. The result is based on the measurement data shown in Fig. 4.2. 


From Data Rate Point of View 

If L-level pulse amplitude modulation (PAM) is considered and there are L equally likely amplitude values placed 

at level ±d , ±‘M . ±(L - 1 )d. In lime reversal system, the waveform filter in the transmitter side is Gr{t) and 

the corresponding frequency response is Gr(f); the time reversal filter is h{-t) and the corresponding frequency 
response is //*(/). So, the average transmitted power is 


cP l? - I r 30 

/ I —-- 

T s 3 


/ OO 

\g t (/)1 2 W (/)l 2 <tf 

■oo 


(4.34) 


In the receiver side, afler matched filter, the minimum distance of the signal is 


- d 1 


‘ /*oo 1 2 

/ \G t (/)| 2 |// (/)| 2 df 

—OO J 


(4.35) 


and the noise power is 


a 2 = K> r pru) i 

^ J - oo 


df 


(4.36) 


where /V 0 /2 is the double sided PSD of the white Gaussian noise. Thus, the PAM probability o symbol error is 


F, 


- 2 ('- 1 Mv 

= 2 0 4 ) q 



3 2T s F t IZo \°T (/)!*!" (/)| 2 df 


P-l N t) \Gr(f)\ 2 df 


(4.37) 


(4.38) 
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where Q(x) is defined as Q(jt) — 
as follows if L is given. 


e 12 dt and the corresponding bit enor rate (BER) H can be determined 



(439) 


Meanwhile we assume that the time duration of Gf{t) is T and T$ - T -f So, the number of bits that can be 
transmitted in each 7 s is 

b — log 2 L (4,40) 


and L can be equal to 2 , 2 s , 2 ; \ — 


In the lime reversal chirp system, Gr(t) is c^(t ) defined in Eq. 4.18. If 7\ f% and L c are given, we can get 27iPi/Wo 
from Eq. 437 and Eq. 4.39 and the data rate of the time reversal chirp system is 



(4.41) 


In the time reversal modal modulation system, there are many parallel orthogonal channels and each channel corrC’ 
sponds to the orthonormal eigenfunction (£) and eigenvalue Kn from Eq. 4.14. We assume kj > > ■ ■. 

Meanwhile, the bits transmitted in these different channels are assumed to be totally independent. 


For the nth channel, the PAM level is I**; Grn(t) is x n 4> n (t) and x n is the value to be determined later; the 
transmitted powei is P tn ; BER is P^; the PAM probability of symbol enor is P n ; the number of bits transmitted in 
each 7 5 is b n , So, from Eq, 434, it follows that. 



and from Eq. 438, it follows that. 
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and 


and 



Pe n 

(4.51) 

log-2 Ln 


tr 

a 

(4.52) 


Thus, ihc total power P t used for time reversal modal modulation is 


Ft = J^P t „ 

n—\ 

d 2 (Ll - 1 ) 

T *k 3 


and the total bits transmitted in each Ts is 


oo 

b = Y, bn 

n= 1 

oo 

= Y1 lo ^ Ln 

n~] 


So, ihc data rate of lime reversal modal modulation system is 



(4.53) 


(4.54) 


(4.55) 


(4.56) 


(4.57) 


From the data rate point of view, we have the optimization problem shown as below. 

(50 

max b — ^2 b n 

n=l 

oo 

S.t. P t =Y 

n=l 

Pbn — constant 

The data rate of time reversal modal modulation system is b/Tg* 


(4.58) 


The iterative loading algorithm shown in Appendix C is used here to solve the above the optimization problem. So, 
the data rale gap is defined as, 




Rm 

R c 

00 

1°S2 I'n 
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(4.59) 


(4.60) 


When BER P b is set to 10“ 4 and the same pulse amplitude modulation (PAM) constellation is considered, the data 
rate gap between time reversal modal modulation system and time reversal chirp system is calculated as follows: 
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When 4 = 2, g rate = 0,3.0,7.0 <dB) for T = 1,10,100 (ns); when L c = t, g rMe = 5.4,7.4,10*2 (dB) for 
T= 1,10,100 (ns); when L c - 8, g ra ie = 7.5,11.1, L3.9 (dB) for T= 1,10,100 (ns). 

From the comparison, we know the performances of lime reversal modal modulation system are better than those 
of time reversal chirp system when Ts > T + 27 V If we want to fully employ the advantage of chirp signal, the 
requirement of Ts should be relaxed. In this case, the introduction of 1S1 is inevitable, in Section 4.3.5, we will 
consider the system performance from BER point of view when Tg is less than 7\ 

4.3*5 BER Performance for Time Reversal Chirp System 

In this section we consider the time reversal chirp system with the existence of ISI. In this case, duration of symbol 
is far less than 1\ Supposing PAM signals are used in a single-user system, the transmitted signal can be expressed 
as 

ao 

s(t) = Y InP ^ - nTs) (4.6!) 


where 




p(t) = Ch(t) * h(-t) 


(4.62) 


l n represents the nth discrete information symbol with duration T$, c^(t) is the chirp signal with the bandwidth W 
and the time duration T, and h(~t) is the time reversed channel impulse response h(t). The received signal is 

OC 


r(t) = Yi l nV(t ~ nT s ) + n(t) 


(4.63) 


{) 


where 


n(t) is AWGN. 


y{t) - p{t) * h(t) 


(4.64) 


The received signal r(t) passes the matched filter followed by a sampler with a sampling rate of 1/7s- The impulse 
response of the matched filter is So the output of the matched filter is expressed as 


oo 

q(t)=Y, r nRpp(t-n^) + v(t) (4.65) 

n=0 

where Rpp(t') is the autocorrelation of p(f) and is the response of the matched filter for AWGN noise n(t). 
Denote 

Rk-n = Rpp{t - nTs) |t=*r s (4.66) 

v k = w(t)[i=*3’ s (4.67) 

wc get 

Qk ~ Inftk-n + Vk (4,68) 

n—U 

= InRk-n + v k 

n—Q,n^k 


(4.69) 
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4 meters 



Figure 4.11: Spectral efficiency in a rectangular metal cavity. 


where I k term represents the expected information symbol of the kth sampling period, the 2 xd term is ISI and v k is 
the additive Gaussian variable at the fcth sampling point. 

Here we performed the time reversal chirp simulation using the channel depicted in Fig, 43 , The time reversal 
chirp system has 2-level PAM and 500Mbps data symbol rate. The frequency range of the chirp signal is from 3GI Iz 
to 10GHz and the duration is 100ns. Perfect synchronization was assumed. The simulation result is depicted in 
Fig, 4 . 10 . From the BER curve we can see that the £fc//V 0 penalty for 10 3 is approximately !7dB, As 1%/Nq 
reaches 25dB, the trend of the slope becomes flat and /i/A o approaches between 10” 4 and 10“ 5 . For a large class 
of UWB applications, it's sufficient for BER at 10 -3 with data rate up to 500MHz in harsh environment with I us 

delay spread. 


4,4 Capacity in Rectangular Metal Cavity 


The channel transfer function is H(f) with bandwidth W = f x - / 0 where f 0 is the starting frequency and /i the end 
frequency. Here, //{/) contains path loss and multipath. The functional relation between path loss and multipath is 
complicated to parameterize. Instead, we use the measured H (/) for calculation. No fading is observed for a fixed 
transceiver pair. 

In the transmitter side, different precoding schemes are employed to shape the spectrum of the signal. Here, water 
filling, time reversal, channel inverse, and constant power spectrum density (PSD) of signal arc considered. Sec 
Appendix A for details. Their relative spectral efficiencies arc compared in Fig. 4J L By simulation it is easy to 
find that the curve of constant PSD can represents that of ideal chirp signal, white spectral efficiency of time reversal 
is almost the same as that of ideal time reversal chirp signal. The channel inverse scheme is good for detection, but 
performs poorly in terms of capacity. 

The capacity is limited by the PSD of the transmission signal. In the US, the FCC requires that the maximum PSD 
be “41,3d Bm/MHz. If the PSD of the AWGN noise is -1 L4dBm/MHz, then the ratio of the transmitted signal 
power and the received noise power, called TX SNR (defined in Appendix A) in Fig. 4,11, is around 70dR. 
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Appendix A Spectral Efficiency 

The PSD of the signal in the transmitter side and noise the receiver side are, respectively, denoted as S(f) and^P, 
The power of the signal is defined as 

P = / 1 S(f)df (4.70) 

Jfn 

The power of noise in the receiver side is expressed as 

* = T (/l " /o, = T w H ' 71) 

with ratio of the transmitted signal power to received noise power (TX SNR) 



(4.72) 


The capacity is 


and the spectral efficiency is 





Nsi 

2 



(4.73) 


(4.74) 


If water filling is used, then 

*11 

S{f) = {t * ~ \H[fW )+ <4J5) 

where (a;) + = max[0, x\, the constant fi is the water level chosen to satisfy the power constraint with equality 



(4.76) 


and the spectral efficiency in this case is 

C /j>g ; 2 (^))+tf 

w fi-fo 


(4.77) 


When lime reversal is used, it follows that 

S{f) = a\H{f)\ 2 

the constant a is the factor chosen to satisfy the poweT constraint with equality 




(4.78) 


(4.79) 

(4.80) 


(4.81) 
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and 


so 


The spectra] efficiency in this case is 
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Kor channel inverse transmission scheme, it follows that 

5(/) = IWjF 

where the constant a is the factor chosen to satisfy the power constraint with equality 
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and the spectral efficiency in this case is 
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(4.8.1) 
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if the transmitted signal is of the constant PSD from j£j to f\. 



so the spectral efficiency is 
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/£ iog,(i + £ffw 

f\ “ /o 
1 1 “ fo 


(4 95) 


( 4 . 96 ) 

( 4 . 97 ) 


Appendix B Maximum SNR after Matched Filter for Time Reversal Modal Modulation 


From Eq, 4.14, we assume the maximum Kn is and the corresponding orthonomia! eigenfunction is 4>(t ), 0 < 
t < i\ 


If the transmitted signal is x<fr (/) * h (—£), then the energy of the transmitted signal is 
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( 4 . 98 ) 

(4.99) 

( 4 . 100 ) 

(4.101) 

(4.102) 

(4.103) 

( 4 . 104 ) 


In the receiver side, after matched filter, the maximum SNR is 


SNR 


max 


[/< 


0 (x4>(t))dt 
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( 4 . 105 ) 

( 4 . 106 ) 

( 4 . 107 ) 

(4.108) 
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Appendix C Iterative Loading Algorithm 


The iterative loading algorithm is employed here to solve the optimization problem. The algorithm can be expressed 

as: 


1 Initialization: bn = 0, n — L 2* 


Tn — 1 1 H — 1 ? 2 , . 


Ptn = 0, n = 1 , 2 ,... Pin = constant, n — 1 , 2 , — 


2 For each channel n, if one more bit is loaded, the incremental power is calculated as follows, 


&Ptn=Ptn(b n +l)-Ptn(bn) 


where P tn (b n ) satisfies 


Pfm (log 2 ^») 


AT) 


Q 


and P tn (b n + 1) satisfies 


Pbu (1 + log 2 hn) = 2 ( 1 - AS) Q 



2 T,P tn (fc n ) 
I N 0 



W,Ptn(bn+ IK, 


/Vo 


(4.109) 


(4.110) 


(4.111) 


3 We get the channel index j according to 

j — erg min AP tn (4.112) 

n 

oo 

If P tn + AP t j < Pt>bj + 1 bj, 2Lj —Lj, P^ + AP £j —* and go to step 2; otherwise, go to step 4. 

71 — 1 

OO 

4 Tlie number of total bits transmitted is b = ^2 b u . 

n^= t 
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Chapter 5 


Chirp UWB Systems 


Chirp UWB—accepted into IEEE 802,15.4a in December 2006—is a promising technology. Our interest mainly lies 
in solving the problems faced with impulse radio—narrowband interference and synchronization. Range extension 
is our goal. 


5,1 Chirp UWB 

5, LI Challenges for IR-UWB 

Though [R-UWB has the advantages of low complexity, low power, etc., it is facing two challenges, limited range 
and Narrow Band Interference (NBI). 

The range for [R-UWB communications is limited by the transmit power. At present, typical range is less than 
30m. In order to extend the range, an IR-UWB system has to increase its transmit power. However, UWB system 
is a power-limited system. Since IR-UWB pulse has ultra-short period while ultra-low duty cycle, increasing the 
transmit power will increase the pulse's Peak-to-Average Power Ratio (PA PR) largely. Such high PA PR makes it 
difficult to implement in hardware. 

The second challenge for UWB is narrow band interference (NBI). Since UWB covers a large bandwidth, it is easy 
to suffer NBI with a Signal to Interference power Ratio (SIR) less than -20dB. NBI suppression algorithms have to 
be used but the complexity of such a system will be very high. 

Timing and synchronization is the classical problem in an impulse radio, even if non-coherent transmitted reference 
and energy detection is used. 


5,1,2 Range Extension 

In this report we would like to propose another UWB system that has the potential to overcome the above challenges. 
It's Chirp UWB, a chirp spread spectrum (CSS) based UWB system. CSS is originated from the RADAR system's 


7i 
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pulse compression theory. Pulse compression is developed to solve the RADAR system's paradox: maximizing the 
detection range while minimizing the resolution of the detected object. Pulse compression requires signals that can 
be simultaneously spread in time and frequency. Spreading in time keeps the transmit power low while accumulates 
high SNR after pulse compression at the receiver. Spreading in frequency ensures high resolution for the detected 
object. Chirp signal inherently has this dual spreading nature. A typical chirp signal has the form of Equation 5.1. 

2ir yf 0 l ± ft 

/□, T CJ ft are center frequency, signal duration and sweep rate, respectively. From Figure 5.1, we can see that chirp 
signal is a linear frequency modulation (LFM) signal. Figure 5.2 is its amplitude spectrum. Us frequency ’sweeps' 
continuously and its duration can be changed independently. Signal with is an up chirp while '-ft' is a down 
chirp. They are matched to each other. Equation 5.2 and Equation 5.3 are the match filter’s impulse response and 
matched output for Equation 5.1, respectively. 


Ch (0 


— cxw { 



T c T c 

- -4- < t < -4 


(5.1) 


h(t) = y/ifiCOB (2 Iff (/uf =F 


ft I 


2 is fit 


< t < — 

“ “ 2 

(5.2) 

T c < t < T c 

(5.3) 


Figure 5.3 illustrates the compressed pulse. After match filtering, a compressed pulse is generated, which has a short 
duration of 2/ B and a peak of y/TB. As a result, a flat envelope pulse with duration T c is compressed after matched 
filtering. Note that the TB product is the processing gain, which increases the SNR after match filtering. 

This feature can solve the challenge of limited range for UWB, because higher SNR in turn compensates the path 
loss costs by distance. SNR gain equals to TB product of a chirp, which can be expressed in Equation 5.4. 

SNR = ^ - TB, in dB (5.4) 

I\Q 

As a result, the receiver sensitivity can be as low as (Ei//V n - TB) in dB. Since the entire chirp signal is spread 
in lime and frequency, an increased signal power can be achieved by increasing 4 or bandwidth, while maintaining 
the same transmit power and PAPR. This is fundamentally different from IK-UWB. Take a CSS system tor example. 
A Chirp UWB pulse has T = Ifts and B = i GHz, and the processing gain is TB = 30 in dB. BPSK modulation 
is used. Under the FCC power spectrum limit and -114dBm noise level, the transmit power is -4 L3dRm and the 
transmit SNR is 72.7dB. The path loss model of LOS channel [1] is depicted in Figure 5.4. With 30dB processing 
gain, the receiver sensitivity is about -20dB f which means that the communication range can be 1000m with BER of 
ID -3 . In this case, when BPSK modulation is used, the data rate is 1 Mbps. 


5.1*3 NBi Suppression 

Chirp based system also has a simple solution for NBI suppression. Chirp signal’s spectrum can be controlled in the 
time domain. If NBI appears at a certain frequency, that frequency can be filtered out by 'gating' in the time domain. 
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Time 


Figure 5J: Waveform of an up chirp 


The block diagram Is dipictcd in Figure 5.5. Figure 5.6 is a 'gated' chirp. Figure 5.7 is its amplitude spectrum. 
Equation 5.5 and Figure 5.5 illustrate the 'gating' process. 


h 9 (t) 
t 



(5.5) 


It can be noticed that a "gated” chirp is summed up with N continuous chirp signals. Each m represents a frequency 
component. A certain frequency can be “filtered” by setting corresponding n\ to 0. If m * — I for all i, the summed 
signal becomes a normal chirp. 

From the amplitude response we can conclude that such 'gated' chirp acts like a notch filter. However, this notch 
filler is easy to implement. It has been reported that this gating process can be done by the surface acoustic wave 
(SAW) device, which is analog (2], Figure 5.8 and Figure 5.9 shows the NBI suppression result. It can be seen 
that when SIR is -20dB, compared with the normal match filter, the NBI is suppressed successfully after the 'gated' 
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Frequency 



Frequency 


Figure 5.2: Amplitude spectrum of chirp 


SAW. It should be noted that this gating process is implemented by a fast adaptive switching network, which is 
digital. Such analog-digital hybrid system makes it easy to control the spectrum. 


5.1.4 Other Advantages 

In addition to the range extension and NBI suppression ability. Chirp LTWB has the advantage of fast synchronization 
because of its analog match filter and high processing gain. Match filter can be implemented by the SAW device in 
the analog domain. SAW is a passive device. A compressed pulse with processing gain will automatically appear 
after match filtering. Timing acquisition becomes a task of peak searching, which is not difficult to implement. 
Also, chirp signal is ’immune’ to frequency offset. A frequency offset A/ will result a timeshifl St in the match 
filter output. Suppose a chirp signal with frequency offset A/ has the form of Equation 5.6. 


s (t) — cos 



(A + A/) t + 





(5.6) 


Equation 5.7 is the output of the match filter. 
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Time 

Figure 5,3: Waveform of a compressed pulse 
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- y/if.t COR (2tt (/o + A/) f) 


sin (rr (A/ + fit) (7’ c - |f|)) 
2rr (A/ + tit ) 


/;</< r c 


(5.7) 


Equation 5.7 shows that the output is still a compressed pulse. However, comparing it with Equation 5.3, there is a 
timeshifi determined by A/. The relationship between A/ and St can be shown in Equation 5.8. 



(5.8) 


Compared with the sweep rate mu, frequency offset is very small, and the influence can be ignored. For example, 
if we add A/ = 1 MHz to a chirp with p = 500 MHz//ts, the timeshifi is only 2 ns. In all. the tasks of timing and 
frequency synchronization can be accomplished together by simple peak searching. 

Since Chirp UWB satisfies the definition of UWB, it has all the merits of traditional UWB. It has the same multipath 
resolution as IR-UWB. Ranging is not a problem for Chirp UWB, since pulse compression with CSS has already 













SNR at Receiver (dB) 


76 


CHAPTER 5 . CHIRP UWB SYSTEMS 



Figure $.4: Path loss model for indoor LOS channel 



Figure 5.5: Block diagram for gating process 


been popular in RADAR applications. Besides, Chirp UWB can be implemented by the analog SAW device, its 
power consumption will be low and it will have long battery life. As a result. Chirp UWB provides the ability to 
enhance current UWB performance. 
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Figure 5.6: Waveform of a ’gated’ chirp 


5.2 Advantages of CSS over DSSS 


CSS is similar lo Direct Sequence Spread Spectrum (DSSS). Instead of spreading the spectrum using PN code 
like DSSS, CSS uses chirp signals. However, CSS is featured with Doppler immunity, fast synchronization and 
controllable spectrum, which are impossible for DSSS. In RADAR applications, DSSS is also a choice for pulse 
compression. But it is well known that BAKER code is very sensitive to Doppler. DSSS with large bandwidth 
and high processing gain requires very high speed Digital-to-Ana log Converter (DAC) and Digital Signal Processor 
(DSP). Both DAC and DSP’s performances arc limited by the state-of-art and they consume a lot of dissipation. 
For synchronization, DSSS needs to do the digital domain correlation to generate a compressed pulse and detect its 
autocorrelation peak. Tliis process cost a lot of hardware resource and increases the system complexity. It is also 
impossible for DSSS to filter out NB1 as easily as CSS, because the pulse for DSSS is fixed. 


5.3 Implementation 

The key device in the Chirp UWB system is SAW. SAW is widely used as bandpass filter, delay line and pulse 
compression Dispersive Delay Line (DDL). It is estimated that the current world-wide production of SAW is well in 
excess of 1000 million devices per year [3]. For CSS purpose, DDL SAW is exclusively used. A DDL based solution 
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Frequency Hz) 


Figure 5.7: Amplitude spectrum of the 'gated' chirp 


was approved by IEEE 802 A 5.4a WPAN task group in March, 2007. In the system, a pair of SAW devices arc used 
as chirp generator and match filter. Each SAW generates an up/down chirp with 80MHz bandwidth and I/*s duration 
[4], It is reported that under 10“ 3 Bit Error Rate (BER), such communication link can support transmission distance 
up to 9800m with IW transmit power [5]. However, the acceptance of IEEE standart only justifies CSS's perfor¬ 
mance for wideband. CSS's UWB application. Chirp UWB, should be of interest. The slate-of-art performances of 
SAW devices can generate a chirp signal with bandwidth over 500MHz [6, 7, 8] or fraction bandwidth over 25% [9], 
satisfying both UWB definitions. In addition, small size, low cost and mass production arc the features of SAW [6], 
which will be perfect for sensor network applications. Combining CSS with UWB, Chirp UWB system will inherit 
the advantages of both CSS and UWB. 


5.4 The History of CSS 


The idea of pulse compression was proposed by different researchers in different countries after World War II 
[10, 11, 12, 13]. In the late fifties, several major labs published their research results on chirp RADAR [14, 15, 16]. 
The first paper applying chirp for communications other than RADAR was published in 1962 by M, Winkler [17], 
where he suggested using up and down chirps to represent binary 'O' and T, respectively. This is the basics 
of a Binary Orthogonal Keyed (BOK) system, utilizing the quasi-orthogonal property of up and down chirps. In 
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Figure 5.8: Normal match filter output with NBI, SIR — -20dB 



1966, Ilata [18] and Gott [19] independently proposed a BOK system for the High Frequency (HF) band. In their 
proposals, chirp signal’s Doppler Effect immunity was utilized. After Gott’s work, Forgan [20] and Edge [21] had 
finished their research on chirp signals and they all reference some Royal Aircraft Establishment and Ministry of 
Defense reports. It seems that their research was supported by the British military and they had shown great interests 
in CSS. The paper of Zaytsev [22] also showed some interests from the Russian military. In the early seventies, Gott 
and his group published several papers on Low Earth Orbit (LEO) satellites and ground-to-air communications in 
the HF band (23, 24, 25, 26, 27]. However, the idea of using SAW device as a chirp generator was first proposed 
by Bush [28] in 1973. Before that, the only practical chirp signal generator was lumped element filters. Interest in 
CSS became less since late seventies. At that rime DSSS and Frequency Hopping Spread Spectrum (FHSS) was the 
more popular choice because of the development of Digital Signal Process (DSP) algorithms. Only a few papers 
mentioned CSS by applying it to DSSS against the Doppler Effect [29,30,31 ]. At that time, SAW device was mainly 
developed for RADAR pulse compression [32, 33, 34], satellite signal processing [35, 36, 37] and channel sounders 
[38,39]. Interest in CSS fell can be explained by the target frequencies and the SAW implementation techniques of 
that time. The excess delay for HF band is in the order of milliseconds, which means that the chirp signal’s duration 
must be no less than that level. Since the velocities in SAW arc thousands of meters per second, the corresponding 
physical size of SAW would be over 1 meter, which is impractical. 

However, when the frequencies are not limited in HF band, SAW devices found its applications. For indoor appli¬ 
cations, the excess delay is in the order of microseconds, which is practical for physical size of SAW. SAW devices 
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Figure 5.9; ’Gated’ match filler output with NBI. SIR = -20dB 


were first applied in the indoor channels by Tsai and Chang in 1994, using BOK scheme [40], After that, several 
groups look interests in CSS. Hucmer [41] in 1998 published another BOK based system. Guglers group [42] and 
Pinkney [43] focused their CSS research on Direct Modulation (DM), another CSS modulation scheme which is 
more efficient than BOK. In their study, the data rate could be as high as 70Mbps with 80MHz bandwidth. In March, 
2007, IEEE approved Nanotron’s CSS physical layer (PITY) in its new wireless standard 802.15.4a. A pair of SAWs 
are used in its BOK scheme. The target applications for such system are Real Time Location Systems (RTLS), in¬ 
dustrial control, sensor networking and medical devices in the ISM band, taking the advantages of CSS; long range, 
high robustness and low power consumption. 


5.5 A Time Reversed Chirp UWB—Forward Looking 


Chirp UWB radio can be regarded as a compromise between the optimum waveform and ease of implementation— 
analog or digital, by selecting chirp waveform for signal modulation. The research proposed here has a significance 
of practical value, and the chirp system is used as a benchmark to be compared against the optimum modal modula¬ 
tion. 

A RF harsh environment shown in Fig. 2.1(a) has been investigated. One simulation example is included here. A 
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Figure 5 JO: The photo of a pair of SAW filters. The size is about l x 2 inches. 


data rale of 500Mbps can be supported using time reversal chirp waveform with a bandwidth 7GHz and waveform 
duration T — 100 ns—the time-band width product is T B — 700. BPSK modulation is used. For a bit error rate 
(BER) of ID -4 , the ratio of bit energy to the noise, E^/Nq, is required to be larger than 23 dB. 

One PhD student (Peng Zhang)—who did his MS thesis on chirp UWB for three years—joined the PFs team in 2006. 
He is currently working on the testbed for chirp UWB* using a SAW chirp device. On this topic, Mr. Zhang is work¬ 
ing with one R&D engineer Dr. N. Guo and another MS student Amanpreet Singh Saini. A pair of SAW filters shown 
in Fig. 5 JO arc available at the Pi’s lab. The vendor is Sawcom Tech, Inc. (http : //www.aawcomtech.comj). 

This is the first phase of the chirp UWB radio. The goal is to prove the concept as quickly as possible. The second 
phase of this project is to cope with the cognitive radio concept, which is the proposed topic for this proposal. The 
SAW filters can be electronically-controlled. Narrowband interference suppression is an important research task, 
due to the co-existence of UWB and narrowband systems such as Wi-Fi and Wi-MAX. 

A chirp UWB system combines the concepts of CSS and UWB with the following features: 

1. Long range transmission 

2. NBI suppression 

3. Fast synchronization 

4. Low transmit power, low probability of intercept 

5. Low power consumption 

6. Low system complexity, low cost 

In summary, chirp UWB is practical and useful. Its combination with time reversal is a promising transmission 
scheme. Time reversal deals with the extremely quasi-continuous multipath—unique to UWB—while a pair of 
chirp pulses is used for signaling. Time-bandwidth-product (TB) solely determines the system performance. The 
range and data rate can be irade-ofTed, by conveniently controlling the time-handwidlh-producl. 
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Chapter 6 


1 VVB Channel Sounding 


In this chapter, UWB channel sounding is discussed. In general, our Lab has the capabilities of sounding the channel 
in both time domain and frequency domain, with and without multiple antennas. The goal of these measurements 
is to evaluate the non-fading wireless communications enabled by time reversal transmission. Both space and time 
are investigated as a function of signal bandwidth. A critical question is: What is the optimum signal bandwidth? 
This question depends on the specific transmission scheme—we are interested in two schemes: time reversal modal 
modulation, and time reversal chirp UWB. 


6.1 Measurement Setup 


The knowledge about the UWB radio channel can be gained by using either time or frequency domain measurement 
technique- We will discuss both techniques in detail in this section* 


6,1.1 Time Domain Measurement Technique 

We will discuss the time domain measurement technique based on impulse transmission. In this technique a narrow 
pulse is sent through the propagation channel and the channel impulse response is recorded by using a Tektronix 
Digital sampling oscilloscope (DSO), The pulse sent through the channel is modulated by a special waveform like 
sine, square or ramp which enables us to analyze the effects of different paths on the received signalp ]* The data 
stored in DSO can be processed by using a computer* The bandwidth of received signal depends on the shape and 
width of transmitted pulse. The time domain measurement setup consists of a pulse generator, transmitter antenna 
and receiver antenna, a triggering signal generator. Low noise amplifier (LNA) and the DSO* The setup for time 
domain measurement is shown in Fig. 6* I 

The whole setup consist of two section, he transmit and receive part. The signal generator and pulse generator 
constitute the transmitter part and DSO along with LNA constitutes the receiver part. The signal generator is used 
to trigger the pulse generator and pulse generator generates the pulse that is transmitted through the channel. On 
the receiver side the signal is passed through LNA and amplified. The final results are displayed on the DSO* 
A triggering signal from signal generator is used to synchronize the DSO to record the measurements. By using 
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Figure 6.1: Time domain measurement setup block diagram 


averaging the signal to noise ratio (SNR) is improved. 

The main advantages of this technique are less complexity, lower cost and channel responses is readily available in 
time domain. With multiple antennas, our system can sound the UWB MIMO channel in the time domain, by using 
two sample modules from Tektronix. 


6.1.2 Frequency Domain Measurement Technique 

In the frequency domain measurement technique, the channel sounding is earned out by sweeping a set of narrow- 
band sinusoid signals (tones) through a wide frequency band. Each tone has the same power (typically 3dBm is used 
in our experiments). The channel frequency response is recorded by using a vector network analyzer (VNA). The 
VNA is operated in transfer function mode where one of its port serves as transmitting port and the other as receiving 
port. S-parameters are used to measure the complex frequency domain transfer function. The setup for frequency 
domain measurements is shown in Fig, 6.2 



Figure 6.2: Frequency domain measurement setup block diagram 


Generally parameter is used where port 1 and 2 are acting as transmitting and receiving port respectively. 
VNA sends a frequency tone / through the channel and the frequency response of the channel is represented as 
821 (f) corresponding to frequency tone /. By sweeping a frequency range from jf) to f\ we can get the frequency 
response of that particular band. IfN is the number of frequency points per sweep with frequency step k Mllz, their 
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relationship with bandwidth B MHz can be represented as 


where B is 



B 

N - 1 


H = fi- fo 


The maximum detectable delay of the channel can be calculated as [2] 


'max — 


N- I 
B 


( 6 . 1 ) 


( 6 . 2 ) 


(63) 


The main advantage of the frequency domain measurement technique is that the cables and connectors can be 
calibrated before measurement to compensate for different kind of errors and frequency dependent losses that can 
occur during the measurement* There are three different kind of measurement errors [3] 


• Systematic errors 

• Random errors 

• Drift errors 


Systematic errors are caused by imperfections in the test equipment and test setup* These errors can be removed by 
calibration* Systematic errors are related to signal leakage, signal reflections, and frequency response. 

Random errors are random in nature and are unpredictable* As a result they cannot be removed by calibrations* 
Random errors are present mainly due to instrument noise* So they can be removed by increasing the source power 
or by narrowing the IF bandwidth. 

Drift errors occurs when the system performance changes after the completion of calibration process. Drift er¬ 
rors occurs mainly due to change in temperature of measurement environment and can be removed by additional 
calibrations* 

A mechanical calibration kit consists of a set of physical devices called standards [4], Each standard has a precisely 
known magnitude and phase response as a function of frequency* These standards are 


• Gpen,short/ofTset shortjoad/sliding load 

• Adapters 

• Precision offsets waveguide or coaxial 


There are various calibration method.The four most widely used are 


• short, open, load, and through (SOLT) 


• LRM 
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* IRRM 

• TRL 


Out of the above mentioned methods the most dominant one is SOLT. With SOLT, the four known references 
are: short, open, load, and through* Once the calibration is done, the calibration state can be saved and recall 
while performing the experiment By using frequency domain measurement technique we can get larger dynamic 
range which improves the measurement precision. The channel impulse response (CIR) which yields the required 
information to characterize the UWB channel, is obtained by taking the inverse fast 1 ourier transform (IFFT) of the 
recorded signal* So this technique requires extra signal processing. 

One of the main limitation of frequency domain measurement technique is that it need static environment through out 
the measurement. But Time domain measurement technique can support non-stationary channels* Another problem 
is the high attenuation of the cables so frequency-domain measurement is good for short range. 


6.2 Measurement Examples 

In this section, we will discuss some of the experiments to gain more insight about the channel sounding* 


6.2,1 Measurement for S1SO Case 

The measurement is performed is for SISO case in three different environments i*e rectangular metal cavity, hallway 
and office to analyze the characteristics of the UWB channel in these environments* The VNA used is Agilent 
N5230A (300kHz-13.5GHz). The calibration method used is SOLT. Table 6.1 list the main parameters lor the 
measurement This set of parameters is used for measurement in all the three cases* I he measurement is done for 


Table 6* 1: Measurement Parameters setup 


Parameter 

value 

Frequency Band 

IGIIz-IOGHz 

Bandwidth 

7GHz 

Number of points 

7001 

Transmission power 

10dBm 

Frequency step 

I MHz 

Antenna Polarization 

vertical 

Averaging number 

128 


Line of sight (LOS) situation. The distance between transmitter antenna and receiver antenna varies from 0.5m to 
4m in steps of 0*5 m and both the antennas are placed at 1 *35 meters above the ground for all the three measurement 
environments. 
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Measurement in Rectangular Metal Cavity 


The rectangular metal cavity used in this measurement measures 16 feet by 8 feet by 8 feet. The material of inside 
walls of the cavity is aluminum. The setup for measurement in rectangular metal cavity is shown in Fig. 6.3 


4 
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Figure 6.3: Setup for measurement in rectangular metal cavity 


The transmitter antenna is fixed, and the receiver antenna is moved from distance 0.5m to 4m in steps of 0.5m w.r.t 
transmitter antenna along the middle line of the rectangular metal cavity. The frequency domain response at distance 
lm, 2m, 3m and 4m between transmitter antenna and receiver antenna is shown in Fig. 6.4 
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Figure 6.4: Frequency domain response //(/) at distance ]m,2m,3m and 4m between transmitter antenna and 
receiver antenna in rectangular metal cavity 


Alter the channel response is measured, zero padding is done between OHz and 3 GHz, conjugate reflection between 
-10 GHz and OHz and inverse fast Fourier transform (IFFT) is used to get channel impulse response in time domain. 
The time domain response at distance lm T 2m, 3m and 4m between transmitter antenna and receiver antenna is 
shown in Fig, 6,5 
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Figure 6.5: Channel impulse response h(t) at distance im, 2m, 3m and 4m between transmitter antenna and receiver 
antenna in rectangular metal cavity 

The channel impulse response forms a channel autocorrelation function according to 

Rhh = h(t ) * h*(-t) (6,4) 

where is the time reversed conjugated version of The autocorrelation of channel pulse response at 

distance I m, 2m, 3m and 4m between transmitter antenna and receiver antenna is shown in Fig. 6.6 
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Figure 6,6: Autocorrelation function Khh(t) at distance Im, 2m, 3m and 4m between transmitter antenna and 
receiver antenna in rectangular metal cavity 


Measurement in Office Environment 


The same procedure of measurement is repeated in office environment. The same set of parameters as shown in 
table 6.1 is used. The frequency domain response at distance i m, 2m, 3m and 4m between transmitter antenna and 
receiver antenna is shown in Fig, 6,7 
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Figure 6.7: Frequency domain response H(f) at distance lm,2m,3m and 4m between transmitter antenna and 
receiver antenna in cilice 


The time domain response at distance Im, 2mu 3m and 4m between transmitter antenna and receiver antenna is 
shown in Fig. 6.8 
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Figure 6.8: Channel impulse response h(t) at distance lm,2m,3m and 4m between transmitter antenna and receiver 
antenna in office 

The autocorrelation of channel pulse response at distance Im, 2m, 3m and 4m between transmitter antenna and 
receiver antenna is shown in Fig. 6.9. 
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Figure 6.9: Autocorrelation function /^(f) at distance lm, 2m, 3m and 4m between transmitter antenna and 
receiver antenna in office 


Measurement in Hallway Environment 


The same procedure of measurement is repeated in hallway environment. The same set of parameters as shown in 
table 6.1 is used. The frequency domain response at distance lm, 2m, 3m and 4m between transmitter antenna and 
receiver antenna is shown in Fig. 6,10 
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Figure 6.10: Frequency domain response H(f) at distance lm,2m,3m and 4m between transmitter antenna and 
receiver antenna in hallway 


The time domain response at distance Im, 2m, 3m and 4m between transmitter antenna and receiver antenna is 
shown in Fig. 6.! I 
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Figure 6.! t: Channel impulse response h(t) at distance Im,2n%3m and 4m between transmitter antenna and receiver 
antenna in hallway 


The autocorrelation of channel pulse response at distance 1m, 2m, 3m and 4m between transmitter antenna and 
receiver antenna is shown in Fig, 6,12 
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Figure 6.12: Autocorrelation function at distance Im,2m,3ni and 4m between transmitter antenna and re- 

ceiver antenna in hallway 

Comparison of Energy for Rectangular Metal Cavity, Office and Hallway Environments 

The energy of channel impulse response for the three eases can be calculated by using 



(6.5) 


where Th is the time duration of the channel. The energy for the three eases is shown in Fig. 6.13 

From the Fig. 6.I3 T it can be seen clearly that the total energy of the channel impulse response in a confined metal 
cavity decays very slowly, as compared to those in olFice and hallway. 
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Figure 6.13: Energy of channel impulse response for rectangular metal cavity^office and hallway environments 


6.2.2 Measurement for MIMO Case 

The measurement is performed for a 3x3 MIMO case using virtual array inside the rectangular metal cavity. The 
measurement parameters used are the same as described in table 6.1. The distance between the antennas in transmit¬ 
ter array and receiver array is 2 feet. The setup for the measurement is shown in Fig. 6.14 
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Figure 6,14: Setup for 3x3 virtual M1MO ease in rectangular metal cavity 
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Figure 6.15: Frequency domain response 11(f) from tansmittcr antenna 1 to different receiver antennas 
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Figure 6J6: Channel impulse response h(t) from transmitter antenna 1 to different receiver antennas 


The frequency domain response and time domain response from transmitter antenna 1 to different receiver antennas 
is shown in Fig. 6.15 and Fig. 6J6 
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Figure 6.17: Frequency domain response //(/) from transmitter antenna 2 to different receiver antennas 



Figure 6.18: Channel impulse response h(t) from transmitter antenna 2 to different receiver antennas 


The frequency domain response and time domain response from transmitter antenna 2 to different receiver antennas 
is shown in Fig. 6J7 and Fig. 6.18 

The frequency domain response and lime domain response from transmitter antenna 3 to different receiver antennas 
is shown in Fig. 6.19 and Fig, 6,20 
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Figure 6.19: Frequency domain response // (/) from transmitter antenna 3 to different receiver antennas 



Figure 6.20: Channel impulse response h(t.) from transmitter antenna 3 to different receiver antennas 


6.3 Conclusion 


In this chapter, channel sounding for UWB channel is discussed. Time domain measurement technique and fre¬ 
quency domain measurement technique are discussed in detail. Then two measurement examples are considered. 
Firstly we compare the characteristics of UWB SISO channel in rectangular metal cavity, in office and in hallway, 
and then the measurement case for a 3x3 UWB MIMO channel in rectangular metal cavity is also considered. 
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Chapter 7 


(Overview of Testbed 


This report documents some of our recent development work with focus on design and implementation aspect. 

The ultimate goal is to achieve range extension by optimizing system design. Specifically, the overall goal is to 
implement a UWB radio test-bed with waveform optimization. The work of design is to provide a feasible solution 
under some theoretical and practical criteria. We proposed a system architecture relying on continuous-lime {or 
waveform level) time reversal technique, though there are a few other types of optimization rules. One of the key 
components in the UWB time reversal system is the programmable waveform generator based on advanced FPGA 
and D/A conversion technologies, A Xillinx Virtex 5 chip, one of the fastest FPGA chips to date, is chosen to 
perform major digital back-end functions at the transmitter. The ultra-wide bandwidth and ultra-high sampling rate 
pose many design and implementation challenges. Major difficulties are associated with the interfaces between the 
FPGA and mixed-signal devices, timing closure issues in FPGA routing, fixed-point implementation with limited 
hardware resources, etc. 

Management wise, in additional to frequently documenting and periodically reporting, we have come up with a 
manageable mechanism including detailed plans and goals for different periods of times, weekly project meetings, 
monthly project reviews, and work assignments on daily basis. 

Currently the system architecture has been stabilized and wc are working toward the next milestone (see milestone 
list in Table t). 

The rest part of this report is organized as follows. 


_Table 1, Milestones _„_ 

Stable Architecture 8/2007 

Time Reversal System with Cable 1/2008 

Time Reversal System over the Air (Some Field i esi) 6/2008 
Fully Functional UWB System with I'ime reversal 1/2009 

System Field Test 6/2009 

Potential Technology Transfer to US Navy 12/2009 
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Chapter 8 


System Design and Architecture 


With the overall goal and constraints from both theoretical and practical perspectives, a completed top-down solution 
has been given. A very broad range of issues need to be considered- Here just name a few: frequency and band, data 
rate, modulation, synchronization, coexistence, commercial availability of devices, etc. Simply speaking, a feasible 
and advanced design must be made. 


8.1 Bandpass Time Reversal 

From implementation point of view, it is too difficult to build a system working from DC to 500 MHz or higher. 
Instead, a pass-band system is more feasible and flexible. To reduce the requirement for sampling rale, a typical 
technique used with digital signal processing is up- and down- conversion. The carrier phase difference between the 
received signal and the output oflocal oscillator is usually unknown. We will show that quadrature conversion can 
take care of this unknow phase difference. 

For convenient description in the following, complex representation is accepted to represent bandpass signals. As¬ 
sume a bandpass sounding signal j s transmitted over the channel, where uj is the carrier frequency and 

0 { \ is an unknown initial phase. All filtering effect introduced over the radio link is taken into account by p(t). Con¬ 
sider a channel with its channel impulse response (OR) h(t) modeled as tap delay line (implying no consideration 
of frequency selective pulse distortion) 

h(t) ~ ^ - n). (8.1) 

t 

Up on receiving the noisy sounding signal 

r(0 = Y, «ip(* ~ Ti)c ,I " {t-ri)+flnl + «(i), (8.2) 

*. 

I 

the receiver’s down-conversion stage outputs a baseband signal 

y(t) - Y a Mt ~ + n'(0> (83) 
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where 61 is an unknown phase introduced in down-converting. There arc a variety of estimation techniques to 
reduce the impact of noise rf(t) and multiple sounding pulse may be necessary. To describe the central idea more 
conveniently, we assume perfect estimation here* Based on this baseband signal y(t) y there can be a few waveform 
for performing time reversal at the transmitter. Consider the following two optional waveforms: 
option I 

h(t) = - Tj). 

i 

option 2 

p(t) = (8.4) 

» 

The option I needs a great deal of effort in the estimation stage, while the option 2 is extremely simple since it 
does not require any parametric estimation. In particular, the option 2 can work directly with “arbitrary" waveform 
generator with need for pulse shaping* 

Let us consider a simple case of one pulse per symbol and use 


p(T 0 - t)e ^ QiP(T 0 - t - r i y (M ‘ +T ’ ,)-e ° +fll+931 


(8.5) 


as the symbol waveform, where To is a constant to guarantee the transmitted signal being causal. The received 
symbol waveform can be expressed as 

dirtjfcpt'Ji, - t - Ti + T fc )e*f w(t+r *“ Tk)_fl " +fllf<?21 

i k 

= - t)e +*+*>) + Yl <Xia k p(Ta - i - tj + r,O^ (Mt+T, ~ T ‘ l) ~ 0n+fll +*1, (8.6) 





which leads to a baseband version after down-conversion: 

]T] afp{To - ty 9 + ^ Q i a fcP(^b -t-Ti + T k y^ n - Tk) * 8 \ 
i i^k 


(8 7) 


where 

0 — ■ 0 q + 0] + 02 “k ^3 (8*8) 

has been applied. In (8.6) and (8*7) the summations and Yli^k correspond to main lobe and side lobes, respec¬ 
tively* For coherent demodulation, there has to be some mechanism to estimate the unknow phase 0 , but this is not 
necessary for noncoherent demodulation* Also, down-conversion is not required for energy detector based scheme. 


8-2 System Architecture 


The transmitter and receiver architectures are shown in ! ig* 1 and 2, where either one can be divided into four major 
functional parts: RF front-end, mixed-signal, digital back-end and controL/intcrface via PC, 

8.3 Link Budget 

The link budget estimation shown below is based on OOK modulation and free space propagation* 
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Transmitter Architecture 



Figure 8.1: Transmitter Architecture. 


Table II. A Link Budget Example. 


Bit rate 

Minimum required SNR per bit 
Center frequency 
Transmitter antenna gain 
Receiver antenna gain 2 
Distance between the two antennas 
Receiver noise figure 
Implementation loss 
Transmitted power to antenna 
Bandwidth (approximately) 

1 m EIRP spectral density over the whole bandwidth 
Path loss 

Received power the antenna 
Link margin 

Proposed minimal receiver sensitivity 


5 Mb/s 

£f>/No. miri - 16 

4 GHz 
2 dBi 
dBi 

40 meters 
7dB 
6dB 

0.51 dBm (1.13 mW) 
500 MHz 
-68.96 dBm/MHz 
76.52 dB 
72.01 dBm 
6.00 dB 
-78.01 dBm 
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Receiver Architecture 



Figure 8.2: Architecture of a reference receiver design. 






























Chapter 9 


RP Front-Ends and Mixed-Signal Stage 


We mainly relay on off-the-shelf products in our test-bed development. Reported in this section include selection of 
components/moduics and some test results. Major components and modules listed in Table III. 

_Table 111. A list of Selected Components/Modules._ 

Local oscillator center frequency 3.5 GHz or 4 GHz 

Modulator (up-convcrter) bandwidth 500 MHz 

Dual-channel D/A sampling rate IGS/s or higher, with LVDS ports 

A/D sampling rale up to 1,5 GS/s 


9.1 Local Oscillator 


PSA3550C produced by Z-Comrminicalion, Inc. is used as the local oscillator (LO) whose frequency range is from 
3545MHz to 3555MHz. Fig. 9.1 shows the picture of evaluation board of PSA3550C. 

PSA3550C includes ADF4106 made by Analogy Device, which is used as the frequency synthesizer. ADF4106 is 
programmable through P5 port on the evaluation board. When we do the lest, the power supply is set to 6Vdc. The 
amplified RF output P2 is connected to Spectrum Analyzer. P3 is connected to the termination. No attenuations are 
used. The onboard 12MHz reference oscillator is used, so PI is not connected to the external reference oscillator 
temporarily. The REF IN frequency is set to 12MHz. The desired pfd reference frequency is set to lOOOKHz. If the 
desired VCO output frequency is set to 3550MHz, the spectrum of the output is shown in Fig, 9.2. 


9.2 Modulator 


TRF3703-15 produced by Texas Instruments is used as the quadrature modulator. I'ig. 9.3 shows the picture ol 
evaluation board of TRF3703-I5, 

The TRF 3703 is a very low-noise direct quadrature modulator, capable of converting complex modulated signals 
from baseband directly up to RF. The TRF3703 is ideal for high performance direct RF modulation from 400MHz 
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Figure 9.1: The picture of evaluation board of PSA3550C used as the local oscillator. 

up to 4GHz. The RF output block consists of a differential to single-end converter and an RF amplifier capable of 
driving a single-end 50fl load without any need of external components [ l ]. 

When we do the test, we connect the 5Vdc supply to headers Wl and W2; we use PS A355OC to supply the LO 
signal at 3555MHz to i I and terminate J2 with 500 to ground; we use a waveform generator to provide the pulse 
waveform whose width is 20ns as the I input signal and set the common mode voltage on the input signals; we 
connect a Spectrum Analyzer to the SMA connector marked RFOUT (J7) and monitor the output Fhc spectrum of 
the output is shown in Fig, 9,4, 


9.3 Digital to Analogy Converter 


MB86064 produced by Fujitsu is used as the Digital to analogy converter (DAC). Fig. 9.5 is the picture of evaluation 
board of MB86064. Fig. 9.6 shows the GUI of control software for MB86064. 

The Fujitsu MB86064 is a dual 14-bit lGSa/s DAC* delivering exceptional dynamic performance. Each high 
performance DAC core is capable of generating multi-stand, multi-carrier communication transmitted signals. DAC 
data is input via two high-speed LVDS ports. These operate in a pseudo double data rate (DDR) mode, with data 
latched on both rising and falling edges of the clock. To simplify system integration the DAC operates from a clock 
running at half the DAC conversion rate [2], A function block diagram is shown in Fig, 9.7. The reasons why wc 
choose MB86064 as DAC in our test-bed can be summarized below. 

From clock point of view, MB86064 requires a RF input clock at half the DAC conversion rate, with sufficient 
spectral purity to not impact the target analog output performance. Because the DAC cores are clocked on both 
rising and falling edges of the input clock, the requirement of the external clock source is reduced. The sampling 
rale of DAC cores is IGsa/s and the frequency of the RF input clock is 500MHz which is easier to be generated by 
the current FPGA devices. At first, we use AD9734 as DAC in our test-bed. Though the sampling rate of AD9734 
is l,2Gsa/s, the DAC clock with frequency of 1.2GHz is needed. 
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Figure 9,2: The spectrum of the PSA3550C output when the desired VCO output frequency is set 3550MH2. 



Figure 9.3: The picture of evaluation board of 1RF3703-I5 used as the quadrature modulator. 


From the number of DAC cores point of view, MB86G64 has two DAC cores which can generate two signals for l/Q 
simultaneously. So, we do not need to buy and integrate two DACs with only one core in the lesl-bed, which will 
greatly reduce the difficulty of the system integration. 

From the interface point of view. The evaluation board of MB86064 is designed to accommodate up to six SMA 
connectors-two for the differential dock input and four for the two differential analog outputs. Boards are supplied 
with transformers on-hoard to perform single ended-toniifferential and diffcrcntial-tCHsingle ended conversions [3]. 
So using MB86064, we can easily make one simple adapter to connect the evaluation board of MB86064 and the 
evaluation board of TRF3703-15. 

From the waveform memory module point of view, MB86064 incorporates a waveform memory module featuring 
two 16k point on-chip waveform memories [4], These allow the DAC cores to be driven with user programmed 
waveforms without the need for external high speed, pattern generators which can increase the flexibility of system 
design. 
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Figure 9,4: The spectrum ofTRF3703-15 output. 



Figure 9.5: The picture of evaluation hoard of mb86064 used as DAC. 


Now, we just test the function of the waveform memory module. According to the user manual, we follow die steps 
to power up and start to use the evaluation board of MB86064. A single-ended clock source is connected to SMA 
connector CLK; the single-ended analog outputs on SMAs IQUTA and IOUTB are connected to the test equipment; 
the jumper links LKJ and LK2 are set as required; the PC USB programming cable is connected to the CONTROL 
I/F header; the power supply cable is connected to the POWER header with correct orientation; the power supplies 
are turned on; the PC control software are installed and operated. 


9.4 MAXI08 A/D 


An SLbit monolithic bipolar A/D converter with sampling rate up to L5Gbps is employed to convert the baseband 
signal (analog signal) into digital domain. A baseband amplifier between the low pass filter following the square law 
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Figure 9.6: The GUI of control software for MB86064. 
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Figure 9.7: The function block diagram of MB86064. 


detector and ADC is necessary for signal conditioning. To facilitate digital interface, the A/D converter features an 
on-chip, selectable 8:16 output demultiplexing that slows the 1.5Gsps data rate to 750 mega-words/second ported 
to two parallel, differential 8-bit, low-voltage (PECL) outputs. A programmable high-frequency clock synthesizer 
is used to generate the sampling dock for the A/D converter The output frequency of the dock synthesizer is 
controlled by FPGA, so variable ADC sampling rate is achieved. The number of bits of the A/D converter is 
dynamically adapted to the conditions of the environment. It can downscale the digital baseband complexity and 
total power assumption, A double-data-rate (DDR) interface with FIFO is implemented in FPGA to connect the 
A/D converter to FPGA, For the design of the high-speed interface between ADC and FPGA, signal integrity has 
become a critical issue. Many signal integrity problems are electromagnetic phenomena in nature and hence related 
to the EMI/EMC, There are two concerns for signal integrity - the timing and the quality of the signal. Signal 
timing mainly depends on the delay caused by the physical length that the signal must propagate. Signal waveform 
distortions can be caused by reflection, cross talk and power/ground noise. An interface board has been carefully 
designed to solve the signal integrity issue. 
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F i gurc 9.8: MAX 108 funct ion d i agram. 


9.5 SPI Control 


In this system. There are many devices need to configure and control, such as LO, AGC, DAC. Each one has its own 
way of initialization and configuration. In order to reduce the system's complexity, SPI Bus is chosen to configure 
all these devices. Also it could be implemented in FPGA,which is easy and efficient. 

SPI standards for Serial Peripheral Interface, It is a simple interface that allows one chip to communicate with one 
or more other chips. It can easily achieve a few Mbps (mega-bits-per-secoods). Basically: (I) it is synchronous. (2) 
It is serial. (3) It is full-duplex. (4) It is not plug-and-play. (5) There is one (and only one) master, and one (or more) 
slaves. In more details: (1) A clock is generated by the master, and one bit of data is transferred each time the clock 
toggles. Because SPI is synchronous and full-duplex, every time the clock toggles, two bits are actually transmitted 
(one in each direction). (2) Data is serialized before being transmitted, so that it fits on a single wire. There are two 
wires for data, one for each direction. (3) The master and slave know beforehand the details of the communication 
(bit order, length of data words exchanged, etc.) (4) The master is always the one who initiates communication. 
Only one slave is active at a time. 

A SPI master can communicate with single slave as well as multiples slaves by connecting most signals in parallel 
and adding SSEL lines. SPI requires 4 wires to be used in between the two chips. As you can see, the wires are called 
SCK, MOSI, MI SO and SSEL, and one of the chip is called the SPI master, while the other the SPI slave. Fig. 9.! 0 
shows A SPI master communicate with one slave and three slaves. 

(1) MOSI - Master device data out. Slave device data in (2) MISO - Master device data in. Slave device data out (3) 
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Figure 9.9: MAX 108 development board. 


SCLK - Clock, Generated by Master Device (4) SSE1 ’Slave device enable, controlled by Master device 

When there are three slaves. The master must have three SSEL lines, and activate only one SSEL line at a time. 
Since all three slaves are connected to the MISQ line, slaves that are not activated must not drive the MISO Imc, 

Let’s assume that the master and slave expect 8-bits data transfers, with MSB transmitted first. Fig. 9.1 1 shows how 
would look a single 8~bits data transfer and its the internal interface. 

(1) The master pulls SSEL down to indicate to the slave that communication is starting (SSEL is active low). (2) The 
master toggles the clock eight times and sends eight data bits on its MOSI line. At the same time it receives eight 
data bits from the slave on the MISO line. (3) The master pulls SSEL up to indicate that the transfer is over. 

If the master had more than one 8-bits data to send/receive, it could keep sending or receiving and rc-asscrt SSEL 
only when it is done. 

In our system, SPI bus is used to implement system control. 

For LO, The evaluation board is supposed to come with a cable for connecting to the printer port oi a PC, The silk 
screen and cable diagram for the evaluation board are shown below as Fig. 9.12. We intend to replace the PC by 
FPGA, Using SPl protocol to write ADF4106 registers through the DB-9 interface. The ADF4106 has a simple 
SPI-compatible serial interface for writing to the device. CLK, DATA, and LE control the data transfer. 

For DAC, Its configuration and control are expected to through Fujitsu’s PC USB Programming cable and accom¬ 
panying PC software, as well as download test vectors to the waveform memory module. The programming cable 
interfaces between a host PC’s USB port and the 4-wire serial interface implemented on the device, The connection 
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Figure 9,10: The block diagram of SPI Bus, 
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Figure 9.11: The block diagram of SPI Transfer. 


is shown in i r ig. 9.13. 

Similarly we also intend to configure the device by SPI protocol. To facilitate this. An appropriate connector should 
be fitted to the target board, providing connections to the signals as: (I) SERIAL IN. (2) SERIAL OUT. (3) SERIAL 
CLK. (4) SERIAL EN. Fig. 9.14 shows this connection. 

Other devices could be configured with SPI protocol similarly. 
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Figure 9,12: The diagram of LO interface. 



Figure 9.13: The diagram of DAC Kit 
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Figure 9. (4: The diagram of DAC interface. 
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Chapter 10 


Digital Back-Ends 


Major function played by the digital back-ends include system control/managemcrit, parameter setting. AGC. syn¬ 
chronization, interfacing with PC, as well as modulation and demodulation. 


10.1 System Control Flow 


10.2 Filter Design 

In order to realize time reversal in the test-bed, wc need to build the time reversal waveform generator in the transmit¬ 
ter side. According to the design of the lesl-bed, the FIR filler built in FPGA will be used to generate the waveform 
of time reversal signaling. Meanwhile, FIR filter will be used in the receiver side to accomplish the functions of 
matched filter, correlation and so on in the digital domain. 


10.2.1 The Typical Methods for Filter Design 

The typical methods for the FIR filter design in FPGA are intellectual property (IP) core based method and user- 
defined coding based method. IP core based method means wc use IP core in FPGA to realize the function of digital 
filter. There arc two kinds of IP cores in FPGA can be employed for this purpose. One kind is filter based IP core 
and the other kind is memory based IP core. The user-defined coding based method means we directly write Verilog 
code to build the FIR filter according to the design specification. 

Multiply accumulate (MAC) FIR filter IP core and distributed arithmetic (DA) FIR filter IP core are the frequently 
used filler based IP cores. The MAC FIR IP core implements a highly configurable, high-performance, and area 
efficient FIR filter. All internal data-paths provide full-precision arithmetic to avoid the possibility of overflow. The 
full-prec ision sum-of-products is presented on the data output port. A set of three handshake control signals provides 
an easy-to-use user interface. The MAC FIR filter IP core uses one or more time-shared multiply accumulate 
functional units to service the N sum-of-product calculations in the filler [1]. The DA FIR filler IP core is also a 
highly parameter!zable, area-efficient high-performance FIR filter. The main difference between the MAC FIR filter 
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and ihe DA FIR filter is the latter employs no explicit multipliers in the design, only look-up tables (LUTs), shift 
registers, and a scaling accumulator. 

In its most obvious and direct form, DA-based computations arc bit-serial in natureserial distributed arithmetic 
i SDA) FIR. The advantage of a distributed arithmetic approach is its efficiency of mechanization. I he basic opera¬ 
tions required are a sequence of table look-ups, additions, subtractions and shifts of the input data sequence [I]. All 
of these functions efficiently map to FPGAs. Input samples are presented to the input parallel-to-scrial shift register 
(PSC) at the input signal sample rate. As the new sample is serialized, the bit-wide output is presented to a bit-serial 
shift register or time-skew buffer (TSB). The TSB stores the input sample history in a bit-serial format and is used in 
forming the required inner-product computation. The nodes in the cascade connection or TSBs are used as address 
inputs to a look-up table. This LUT stores all possible partial products over the filter coefficient space. 

In a conventional MAC FIR filter realization, the sample throughput is coupled to the filter length. As the filter length 
is increased, the system sample rate is proportionately decreased. With a DA architecture, the system sample rate 
is related to the bit precision of the input data samples. The trade off introduced here is one of silicon area (FPGA 
logic resources) for lime. As the filter length is increased in a DA FIR filter, more logic resources are consumed, 
but throughput is maintained [!]. Meanwhile, one approach to increase the speed of FIR filter using DA architecture 
is to partition the input words into M subwords and process these subwords in parallel [1]. This method requires 
M-times as many memory look-up tables and so comes at a cost of increased storage requirements. Maximum speed 
is achieved by factoring the input variables into single-bit subwords. The resulting structure is a fully parallel DA 
(PDA) FIR filter. With this factoring a new output sample is computed on each clock cycle. Another advantage of 
the DA FIR filter IP core is the on-line coefficient reload which provides additional flexibility to the system design. 
While the new coefficient values are being loaded, and some internal data structures are subsequently initialized, 
the filter ceases to process input samples, i he coefficient reload time is a function of the filter length and type [1 ], 
When we use the DA FIR filter IP core to build the single rate FIR filter, the reload delay is 210 clock cycles if the 
number of taps is 10, the coefficient width is 10 and the input width is 1. Fig. 10.3 shows the simulation result of 
the coefficient timing using DA FIR filter IP core in FPGA. 

If we want to use memory based IP core to realize the FIR filter, FIFO Generator IP core and Block Memory 
Generator IP core can be employed. The basic thought of this method is that the data of the candidate waveform arc 
put into the memory firstly and then the data are read from the memory according to some control signals later. IP 
core in FPGA has powerful functions. Even though IP core can work efficiently, it will occupy more resources than 
we expect. So if our design is not complex, we can directly write Verilog code to build FIR filter. For example, we 
can build the shift register and adder by coding to generate the waveform of time reversal signaling. 


10.2.2 Quantization Consideration 

Quantization issue is very important for the FIR filter design in FPGA, because any number in FPGA should be finite 
word length. Quantization will introduce error in the filter implementation. There are three sources of quantization 
error - coefficient quantization, input quantization and the internal rounding that is the quantization of the results 
of arithmetic operations within the filter [2]. According to [2], if each coefficient is quantized to m + 1 bits and 
the input is also quantized to n + 1 bits where n < m, the output quantization noise power is nearly optimal for 
internal rounding precision set to less than m + 1 and greater than n + I bits if there is no restriction on the output 
quantization. However, in our test-bed, the output quantization is fixed and set to 14 bits, because the input resolution 
of DAC is 14 bits. Meanwhile if we use shift register to build time reversal waveform generator, we will only have 
coefficient quantization and internal rounding problems. If the lime reversal waveform has no overlapping in the 
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transmitter side, the coefficient quantization should be 14 bits which is the same as the output quantization. Based 
on Dr, Guo’s research results, nearly half of the ideal energy of the channel impulse response can be attained in the 
receiver side using 1 bit coefficient if inter-symbol interference (ISI) is not present Because mono-bit resolution will 
cause the big side lobe in the received signal. So, if waveform overlapping is required in the transmitter side and ISI 
is present in the receiver side, the resolution of the coefficient and the internal rounding should be increased in order 
to suppress the side lobe and increase the energy of the main lobe in the received signal. Because the channel state 
varies in different scenarios, we should make the tradeoffs between the system performance and system complexity 
to determine the resolution of the coefficient and the internal rounding when the output resolution is fixed and given* 
This research issue will be further studied. 

In sum, no matter what method we use, we should follow some steps to design the FIR filter in FPGA. (1) Use 
Matlab to design a floating-point FIR filler, (2) Convert the floating-point FIR filter into fixed-point FIR filter using 
Maltab. (3) Build a hardware description in FPGA and verify that the fixed-point design functions correctly. (4) 
Evaluate the fixed-point FIR filter in terms of its performance and hardware efficiency [3]. (5) Program FPGA and 
verify the FIR filler’s performance experimentally on the workbench. 
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Figure 10.1: System control (low for transmitter. 
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Figure 10.2: System control flow for receiver. 
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Figure 103: The simulation result of coefficient reload timing using DA FIR filter IP core in FPCJA. 
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Chapter 11 


Graphic User Interface (GUI) between FPGA 
and PC 

In this section wc’H introduce the PC interface in our test-bed system. The PC interface is responsible for in¬ 
put/output communications between host PC and FPGA. 


11.1 Architecture 

Fig. 11.1 shows the architecture of the PC/FPGA connection with PC interface. Here wc use laptop as the host PC, 
since it’s powerful enough and easy to carry. The data acquisition (DAQ) board acts as the PC interlace. For digital 
data output, DAQ gets data from PC then send it to FPGA. For digital data input, DAQ samples data from FPGA 
then sends it to PC. We use USB cable for the interface between PC and DAQ while (lying lead wire for DAQ/FPGA 
connection (Fig. 11.2). 


11.2 Communications between FPGA and host PC 


FPGA is extensively used in our test-bed system. However, if we want to change the parameters in FPGA for testing, 
there will be some difficulties. We can let the FPGA change the parameters by itself, but this will take up the valuable 
FPGA resource. Also, FPGA does not have graphical user interface (GUI). Wc must change the codes and download 
them to FPGA each lime, which is time consuming and not intuitive. With a PC interface between host PC and 
FPGA, wc can bypass these difficulties. We can set up the parameters in GUI then transmit them to FPGA via PC 



Figure 11.1: The architecture of the PC/FPGA connection with PC interface 
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Figure 11.2: The flying lead wire interface. 



Figure 11.3: Byte Paradigm, GP22050 DAQ board 


interface without making any changes in FPGA, This is much more efficient and intuitive. 

We can also use PC interface to read data from FPGA. Though this function can be done by logic analyzer, it can 
only display static signals, not real time signals. Also, it does not support real time signal processing. With PC 
interface bridged between FPGA and PC, real time data can be sampled and stored to PC. After capturing, the data 

can be processed and viewed in the monitor in real time. 


11.3 Implementation 


We choose Byte Paradigm’s GP22050 as our DAQ (Fig. 11.3). It is a powerful DAQ supporting up to 16 bidirectional 
data and 6 bidirectional control signals. The maximum processing rate is SOMIIz. Therefore, the maximum buret 
throughput can be up to lOOMByte/s. GP22050 is also provided with C language interface. The C functions can be 
called to control DAQ by Visual C (VC), Labview and other programs [1,2]. Thus, a user-friendly GUI can be made 

to control the PC interface between DAQ and FPGA. 
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Figure 11.4: Data stream structure for PC/DAQ connection. 
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Figure 11.5: The relationship of the clock, trigger and data signals. 


H.4 Application: BER Measurement 


Now we are developing the bit error rate (BER) measurement application. To get the BER, we need to compare the 
received data with transmitted data. To simplify the process, the transmitted data periodic and known at the receiver. 
Also, at the beginning of each period, there is a trigger signal to inform the program when to start comparing. In 
our PC interface based system, the received data, trigger signal and clock signal are the output ot FPGA. DAQ uses 
three data streams to get those signals. The data stream structure is illustrated in Fig. 11.4. DAQ will not start 
sampling until it receives the trigger signal. After 5 clock cycles of the trigger signal [3], DAQ samples data from 
FPGA at each clock rising edge. The relationship of the clock, trigger and data signals is illustrated in Fig. 11.5. 
DAQ samples at rising edge of a clock and it only starts sampling after trigger signal is delected [3], The maximum 
sampling rate is the same as DAQ’s maximum sampling rate, which is SOM Hz [4]. Since DAQ can sample 16 data 
lines in parallel, the maximum throughput can be lOOMbyte/s [4], 

Fig. 11.6 is the BER measurement software GUI. Before BER measurement, we need to initialize the measurement 
environment in the “Settings” dialog first (Fig. 11.7). If “Stop Manually” is chosen, the software will keep running 
until the “Stop” button in Fig. 11.6 is clicked. Otherwise, it stops when it reaches “Number of bits” or “Errors” set 
by user. BER values will be updated every 100ms in the “BER” edit box. A new BER value will be plotted onto the 
"Real time display” box every 1 second. BER VS Distance (or SNR) curves will be displayed on “Statistic display” 

box for each environment. They can be saved and loaded any lime. 

This software is stilt under construction. It is expected to be finished in one week. 
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Figure 11,6: The BER software GUI. 
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Figure 11.7: The Settings dialog. This dialog is used to initialize the measurement environment. 
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11.5 Application: FPGA Test 

Wc arc going to develop an FPGA test software. In order to test FPGA codes, we need to change the parameters 
frequently, such as FIR filter’s coefficients, system data rate, transmitted data, etc.. Since FPGA does not have GUI, 
we can either change the code each time or let FPGA change the parameters by itself. Both of these methods are lime 
consuming and not intuitive. With DAQ, however, FPGA will be able to read parameters from PC's GUI without 
changing the codes. This is a much more intuitive way and saves time. 

GP22050’s “arbitrary digital waveform generator” (ADWG) function [4] can be applied to our FPGA test applica¬ 
tion. The ADWG function is able to generate arbitrary 16 bit-depth output at a rale of 50MHz. As a result, GP22050 
supports setting 16 coefficients simultaneously in FPGA at a high speed. With the ADWG Library [2] provided by 
the device, wc can customize the GUI design and make the test much more intuitive. 


11.6 Limitations 

Although PC interface provides a way for communications between FPGA and PC, the processing speed is limited. 
Since PC utilizes a lot of resources to support the operation system and other necessary softwares, the signal pro¬ 
cessing speed is far lower than the speed of the DAQ. According to result of a test with a Pentium 4 3 GHz PC', the 
maximum continuous throughput is only llMByte/s, which is far less than the burst throughput of lOOMBytc/s. 
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